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ABSTRACT 

This final report under Contract NAS1-13290 documents work done in support 
of a NASA program designed to evaluate the utility of the OMEGA navigation sys- 
tem. With support from the Research Triangle Institute, NASA personnel at the 
Langley Research Center began an extensive experimental program in 1973. Par- 
ticular emphasis was placed on collecting statistically significant quantities 
of OMEGA data to determine through analysis the achievable accuracies with the 
differential and composite modes of operation and to define the major contrib- 
uting factors to differential OMEGA error. One objective of this study has 
been to investigate the potential application of OMEGA navigation to civil 
aviation use. This report is the conclusion of the analysis from which pre- 
liminary results were presented in an earlier final report prepared under 
Contract NAS1-12043. 

Major attention in this report is given to an analysis of receiver re- 
peatability in measuring OMEGA phase data. In this work, repeatability is 
defined as the ability of two like receivers which are co-located to achieve 
the same LOP phase readings. In the context of differential OMEGA this limits 
the achievable position error. In an early phase of the experimental program, 
modal interference of the North Dakota OMEGA signal caused extremely large 
differential errors, particularly during nighttime hours. Specific data anal- 
ysis is presented in this report in support of this conclusion. A propagation 
model is described which has been used in the analysis of propagation anomalies. 
Composite OMEGA analysis is presented in terms of carrier phase correlation 
analysis and the determination of carrier phase weighting coefficients for 
minimizing composite phase variation. Differential OMEGA error analysis is 
presented for receiver separations of up to 600 n.mi. Three frequency analysis 
includes LOP error and position error based on three and four OMEGA transmis- 
sions. Finally, results of phase-amplitude correlation studies are presented. 
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CHAPTER 1 


INTRODUCTION 

Since 1972 the NASA Langley Research Center at Hampton, Virginia 
has been involved with an investigation of the OMEGA navigation system. 
This investigation has been designed to determine applicability of the 
system for general aviation use. An experimental program involving 
collection of OMEGA data and subsequent analysis was initiated in 
September 1973. OMEGA data were recorded according to this plan 
through 1975. The Research Triangle Institute has supported this NASA 
program under a previous contract, NAS1-12043 (ref. 1), and under this 
current contract, NAS1-13290. Associated work involving development 
and analysis of navigation accuracies using a feasibility model of a 
low-cost general aviation airborne OMEGA processor is currently supported 
by RTI under contract NAS1-14005 (ref. 2). 

This report serves to provide a comprehensive analysis of OMEGA 
position estimate errors considering various modes of operation using 
data recorded by NASA personnel. Analysis of receiver system errors 
are also included. Propagation anomalies, including particular emphasis 
on inaccuracies resulting from time-varying modal interference in the 
OMEGA VLF propagation, has been a major concern of this study. This has 
led to development by RTI personnel of a VLF propagation prediciton 
model which is capable of calculating PPC (propagation prediction cor- 
rections) for OMEGA phase considering all significant modes of propaga- 
tion. 

Considerable effort under this contract has been associated with 
supporting NASA personnel in reducing OMEGA data, preparing plots of 
data for visual editing, editing data, and summarizing the data in a 
form suitable for position estimation and error evaluation. Data plots 
were so voluminous that a limited number have been prepared for NASA- 
Langley and are not presented in this report. Much of the LOP analysis, 
differential error analysis, and position error analysis results are 
also not included. These have been summarized in this report in support 



of the conclusions made. 

Use has been made of regression analysis and analysis of variance 
procedures to present results of a parametric analysis of the differen- 
tial OMEGA navigation mode. An attempt has been made to generalize 
these results beyond the particular geographical region in which data 
were accumulated. In completing this analysis, a significant, unex- 
plained error contribution was observed in the receiver/recording sys- 
tem as evidenced by the side-by-side test results. In the experimental 
plan, these were originally intended to allow for possible needed cali- 
bration of the receivers. Ultimately, these showed that random varia- 
tions in receiver performance were a limiting factor in obtaining de- 
finitive results from regression analysis. An understanding of the 
fundamental source of the side-by-side phase difference variations is 
not complete. One hypothesis involving receiver phase measurement 
variation with signal amplitude is discussed. Although the two re- 
ceivers used to record the OMEGA data do not have the same characteris- 
tics, there appears to be little, if any, correlation between variation 
in amplitude differences at the two receivers and measured phase varia- 
tions with the receivers co-located. 

In this report Chapter 2 describes the overall data gathering ex- 
perimental program and analyzes the results of side-by-side tests with 
particular emphasis on receiver repeatability. Also, phase-amplitude 
correlation results are discussed. Chapter 3 discusses the results of 
differential OMEGA studies. These two chapters relate to tasks 1 and 
2 of the contract. In Chapter 4, the VLF propagation model development 
is discussed which relates to task 4 of the contract. Chapter 5 des- 
cribes results of composite OMEGA studies as per contract tasks 1, 2 
and 3. Using actual data and the propagation model, propagation anoma- 
lies are discussed in Chapter 6 relating to contract task 5. Tasks 6 
and 7 of the contract are described briefly in Chapter 2. Under task 
8, two publications which relate to these studies (refs. 3 and 4) des- 
cribe parametric analysis of differential OMEGA and time varying modal 
interference. 

Several appendices are included giving supporting data and analy- 
sis results described within the report. 
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CHAPTER 2 


DIFFERENTIAL OMEGA EXPERIMENTAL PROGRAM 

To evaluate the OMEGA navigation system a two-year data gathering 
experimental program was undertaken. Beginning in September, 1973, two 
OMEGA receiver sets with associated digital recording equipment were used 
to collect phase and amplitude data at all three OMEGA frequencies. With 
a four station capability (A-Norway, B-Trinidad (segment G through March 
1978), C-Hawaii, D-N. Dakota) one of the receivers was maintained in a 
laboratory environment at the Langley Research Center (LRC) , Hampton, 

Virginia. The other receiver system was mounted in a travel trailer and 
periodically moved from Hampton to each of twelve different remote locations 
at varying distances from Hampton. Additionally, six locations extending 
into Florida were used on a one-time basis. Figure 2-1 illustrates the 
various locations used. Several visits to each of the twelve primary loca- 
tions were made during the course of the data gathering. During each visit, 
both the "base" receiver at Hampton and the "remote" receiver were operated 
simultaneously for a period of three to five days. 

This plan was designed for the primary purpose of obtaining data 
suitable for analysis of the differential mode. Of course there are suf- 
ficient data to evaluate other modes including ordinary OMEGA, composite 
OMEGA, and to investigate the effects of modal interference and other 
anomalous behavior. A previous report (ref. 1) has described the plan in 
more detail and also elaborates on the receivers used including an analysis 
of receiver performance. The plan was not complete at the time of the earlier 
report, therefore, this description will include the entire plan. 

In the course of experiment design, two primary objectives were de- 
fined. One was to make measurements to adequately determine the variation 
of differential OMEGA with respect to receiver pair separation range and 
relative orientation. The second was to collect adequate data to gain in- 
sight about temporal variations, i.e., hour-to-hour variation, day-to-day 
repeatability, and seasonal variations. Previous differential tests (refs. 

5, 6 and 7) had examined geometric effects in other local areas, but generally 
the conclusions were limited because relatively few receiver locations were used. 
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Figure 2-2 illustrates the thirteen primary locations used (see 
Figure 2-5 in ref. 1) with the separation range and azimuth from Hampton. 
Table 2-1 is a schedule of site visits for the remote receiver system. 

As described in a previous report (ref. 1), each receiver system includes a 
Tracor 599R OMEGA receiver. The numeral within each darkened area of 
Figure 2-3 indicates the frequency selected on the Tracor receiver (1- 
.10.2 kHz, 3-13.6 kHz). Figure 2-3 displays the distribution of receiver 
separation ranges, azimuths, and the site revisitation period. 

During data gathering the N. Dakota transmission at full power and 
the Trinidad transmission at under 1 kw were the only signals continuously 
available. Norway data were available during about seventy-five percent 
of the time and Hawaii was on the air for at least one visit to each of the 
receiver sites during the latter stages of the experiment. The major 
problems with the OMEGA signals were fairly regular occurrences of modal 
interference on the North Dakota transmission, particularly at night> 

(ref. 3) and poor signal strength on the Norway transmission. 

There were 83 data sets generated through repositioning of the remote 
receiver. Of these, six were not analyzed because of severe problems with 
the recording system or other equipment related problems. Six of the 
remaining sets were for separation ranges which extended beyond the con- 
ventional "differential region." Overall, seventeen sets (22% of the good 
data sets) provide measurements with the receivers co-located or side-by-side 
(SxS). These were interspersed throughout the experiment to allow some 
measure of inherent receiver repeatability (differential repeatability). A 
significant portion of the subsequent analysis has been relative to the SxS 
situation. Some significant random variations in the mean differential 
error between data periods were observed. A detailed analysis of these 
observations follows. 

2.1 Receiver Repeatability 

Receiver repeatability has been used in a strict sense (refs. 8 
and 9 ) to mean capability of a receiver to indicate the same phase measure- 
ments given the same input signal phase at different times . With respect 
to this analysis, receiver repeatability is used in a differential sense. 

The capability of a differential pair of receivers when co-located to 
yield a differential phase measurement which is not dependent upon 
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Hampton, Va. 

0 

0 

— 

YKT* 

Yorktown, Va. 

8.3 

15.37 

328.6 

FEU 

Fort Eustis, Va. 

11.2 

20.74 

286.4 

FIS 

Fisherman Island, Va. 

19.2 

35.56 

90.8 

PGO 

Prince George, Va. 

40.0 

74.08 

277.4 

AHO 

Ahoskie, N.C. 

56.9 

105.38 

213.2 

WAL 

Wallops Island, Va. 

66.0 

122.23 

40.8 

BOD 

Bodie Island, N.C. 

85.2 

157.79 

152.3 

NCC 

Naval Comm Center, Md. 

101.3 

187.61 

347.4 

FOU 

Fountain, N.C. 

103.1 

190.86 

216.3 

RTI 

Research Triangle, N.C. 

139.8 

258.91 

239.9 

BDF 

Bedford, Va. 

151.8 

281.13 

280.5 

THV 

Thayerville, Md. 

200.0 

370.4 

317.4 

CHA 

Chadburn, N.C. 

204.0 

377.81 

216.2 

WAN 

Wando , S . C . 

301.1 

577.64 

215.4 

RMT 

Rich Mountain, N.C. 

326.9 

605.40 

251.26 

MI 

Midway , Ga . 

405.9 

751.73 

219.6 

BLC 

Black Creek, Fla. 

499.2 

924.52 

213.8 

WEE 

Weeki-Wachee , Fla. 

598.9 

1109.16 

212.9 


Figure 2-2. Experimental Plan Receiver Site Locations. 
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Figure 2-3, For Twelve Primary Sites of Remote Receiver: (a) Histogram of Azimuths from Hampton, Va. to 
Various Sites, (b) Histogram of Differential Ranges, (c) Histogram of Site Revisit Delay, 
(d) Cumulative Distribution of Revisit Delay Time. 







time._ When used in this sense, receiver repeatability is determined by 
the capability of each individual receiver to independently indicate the 
same measure of the same phase at different times as well as any differences 
in S curve type error which may exist between receivers. Repeatability 
error is typically random so that in a differential sense both correlated 
and uncorrelated (between receivers) contributions to error are included. 

Any mutual coupling between antenna configurations can also affect repeat- 
ability. 

Receivers used in this experimental plan were co-located (SxS set-up) 
at frequent intervals so that a good estimate of differential receiver 
repeatability could be obtained. Any drift in phase measurement could 
conceivably be accounted for through this ’Calibration” procedure. 

As data have been analyzed, a random repeatability error has been 
observed. This error for any given receiver SxS set-up generally has a 
small variation about the mean (< 1.0 cec) however, there is a rather large 
variation in the mean between periods (^2 cec). The variation in the 
period means as well as the range of the mean values does change with 
frequency. Following is a discussion of this analysis. 

The SxS data have been analyzed primarily on the basis of one-hour 
averages of 10-second sample values of phase. This corresponds to low-pass 
filtering on data recorded from a receiver which has a phase measurement 
time constant of approximately one minute. The resolution of the phase 
measurement circuitry in each digital receiver is 0.25 cec at 10.2 kHz, 

0.28 cec at 11.3 kHz, and 0.33 cec at 13.6 kHz. Figures 2-4, 2-5, and 
2-6 are plots of hourly average differential LOP phase for approximately 
75% of the SxS situations. These plots show full 24-hour period data with 
data sets overlaid. Daytime periods at Hampton, Va., correspond approxi- 
mately to GMT hours 12-20 with exact times dependent on time of year. 

Figure 2-4 shows LOP B-D (TRI-N.DK) differential measurements at all three 
frequencies. More overall spread is evident in the 13.6 kHz data with a 
slightly positive overall bias in the 11.3 kHz data and slightly negative 
at 13.6 kHz. The 10.2 kHz data is rather evenly spread between +1.0 cec. 
Although it is not illustrated very well in these figures, variations in 
day-to-day repeatability are small compared to the overall variations. At 
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Figure 2-5. Continued. 
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Figure 2-6. Continued. 







10.2 and 13.6 kHz there is one transition period during which the differen- 
tial phase error is definitely outside the envelope. This is not explained. 
From these data the 11.3 kHz signal appears to be the most repeatable. In 
Figures 2-5 and 2-6, LOPs A-B (NOR-TRI) and A-D (NOR-NDK) are shown for the 
same data. These show somewhat more overall variation and can be attrib- 
uted to the receiver performance tracking the Norway signal which is 
generally 40-50 dB below N. Dakota at the Hampton, Va. location. These 
low SNR conditions also resulted in typically larger hourly standard 
deviation of differential LOP phase values for those LOPs shown. Any hours 
with phase standard deviation greater than 4.0 cec have been deleted from 
these presentations. No one frequency is consistently best, however, the 
A-B 13.6 kHz nighttime repeatability is particularly good. Figure 2-7 
provides a plot of the Tracor 599R receiver differential data to illustrate 
the relative quality of the NASA digital receiver performance. 

No apparent trends are identifiable in the data set to data set 
variations. This period-to-period randomness is further illustrated in 
Table 2-2 with data set means and standard deviations. Here XMN refers 
to the average hourly mean value for the data set and STD is the standard 
deviation of hourly mean values. 

Data from station C (HAW) was not available until approximately 
December 1974 and has not been used in this discussion of period-to-period 
randomness versus within period repeatability. 

In analyzing repeatability error, of particular concern in the differen- 
tial mode is the effect of this error on position estimation. These same 
data have been analyzed in combination form in terms of relative position 
error. This represents a baseline error contribution which is dependent on 
the receiver system itself and is independent of those parameters which may 
be of greater interest in evaluating the concept: receiver separation 

distance and azimuth; phase dispersion; and local conditions. Position 
error here is used in the sense of east-west (x) and north-south (y) 
variation from the true position assuming that the differential LOP phase 
errors are used to estimate a position fix relative to the true position. 

Let 


A-f - [A<t> x , A<1> 2 , . . - A({H 
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Table 2-2. Side-by-Side Average Hourly Mean and Hourly Mean 
Standard Deviation by Data Set (Full Period) . 



10.2 


XMN 


NOR-HAW 


11.3 


13.6 


STD 


0.53 2.25 1.34 2.78 0.85 2.16 


3.13 6.43 1.54 
0.72 1.96 3.41 

0.12 I.43I 3.42 
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4l|5.05h .38 1.78 
42i 5.02J- .44 1.47 


.12 1.43) 3.42j 5.02|- .44 1.47 
.33 1-52)- .04| 0.97: 0.48 1.37 
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0.54jl.06 0.77 
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0.59 
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2.31 
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2.37; 
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4.36 
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0.44 

1.63 

0.22 j 1. 27: 3.34 

: 4.91 

0.60 
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0.95 

2.10 
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30j .77 
,51) 0.57 
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0 • 16 ! 0 . 8l| 0.85 
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represent a vector of differential LOP phase measurements from the SxS 
data. These can be described in terms of a vector e = [e^, e^] = [x,y] as 


where 


A<p = [M] e 


[M] = 


Y 1 cos8 1 -Y 1 sin8 1 
Y 2 cos8 2 -Y 2 sinB 2 


Y cos8 -Y' sing 
n n n n 


(2-1) 


with 8^ the CW angle which is 90° less than the angle between north and the 
LOP gradient at the true position. The y ^ values are LOP gradients in cec/km 
at the respective frequency of consideration. Considering stations NOR, TRI 
HAW, and NDK, the 8^ and y^ values at Hampton, Virginia, are given in 
Table 2-3. 

Given (2-1) the least squares estimate of position error is 

7 = tM T M] _1 M T A7 . (2-2) 

Since only LOP data are available, the vector A<j> can never have more than 

T -1 

three linearly independent elements and [MM] is 2x2. Table 2-4 summarizes 
(2-2) for all the combinations of three linearly independent LOPs taken from 
four transmitter phase measurements at the Hampton site. 

For the situation when only three transmitter phase measurements are 
available, there can be only two independent LOP measurements, [M] is 2x2, 
and (2-2) reduces to 


s' = [M] -1 A<j> 


(2-3) 
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Table 2-3. LOP Angle and Gradient Values at Hampton, Va. 


LOP 

g (degrees cn) 

Y 

10.2 kHz 

Gradient Values (cec/km) 
Y 

11.3 kHz 

Y 

13.6 kHz 

* 





AB 

89.07° 

5.86 

6.50 

7.81 

AC 

155.76° 

5.43 

6.03 

7.26 

AD 

167.70° 

4.44 

4.98 

5.97 

* 

BC 

216.09° 

6.25 

6.94 

8.35 

* 

BD 

227.93° 

6.68 

7.43 

8.92 

CD 

294.62° 

1.38 

1.57 

1.86 


*B is TRIN 
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Table 2-4. Position Estimate Trans forniat ions at 'Hampton, Va. 
for 10.2 kHz, 4 Transmitters 


Position 

Error 

AB* 

Differential 

AC 

LOP Phase 
AD 

Measurement. Error 
*BC *BD 

CD 

E 

X 

-.0229 

-.1266 


-.0845 


E 

y 

-.0960 

-.0576 


.0618 


E 

X 

.0778 

-.1995 



.0110 

E 

y 

-.1636 

.0001 



.0278 

E 

X 

-.0424 

- 

.1165 

•.0985 


E 

y 

-.1329 

- 

.0547 

.0448 


E 

X 

.0415 


.2254 


.0231 

E 

y 

-.1629 


.0006 


.0315 

E 

X 

-.1083 



-.1941 

.0553 

E 

y 

-.1583 



.0023 

.0459 

E 

X 

-.1636 



-.2170 

.0844 

E 

y 

-.1557 



.0034 

.0513 

E 

X 


-.1187 

- 

.0795 

.0211 

E 

y 


-.1510 


.1563 

.0707 

E 

X 


-.1378 


-.0684 

.0244 

E 

y 


-.1174 


.1367 

.0492 

E 

X 


- 

.1557 

.0580 

.0502 

E 

y 


- 

.1920 

* 1802 

.1224 

E 

X 


- 

.1753 

-.0465 

.0492 

E 

y 


— 

.1467 

.1538 

.0785 


* B is Trinidad 
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Table 2-5 summarizes the transformations applicable to (2-3). Given 
a set of differential LOP measurements which represent all of the possible 
combinations of LOP values from a given transmitter set, any of several 
transformations can be used to get a position error estimate. The position 
estimate will be independent of which LOP pair is selected. Station 
triplets are selected on the basis of ability to track phase and geometric 
coverage. 

Figures 2-8 through 2-13 provide three and four station SxS hourly 
position error plots based on the measured differential LOP errors. Also 
shown on these figures is the overall RMS error for the hourly values and 
the CEP (circle of equal probability) radius for each set. Data for all 
three OMEGA frequencies is illustrated. It can be noted that the most 
significant spread in the position errors is generally in the direction of 
station A (NOR). The signal received from A has the largest phase variation 
and no other available station provides good geometric coverage in the 
direction of A. 

In analyzing these results, the differential RMS position error is 
consistently within 0.5 km. Generally the 13.6 kHz derived position error 
is the best. The maximum resolution error is on the order of 0.1 km and 
is dependent upon the transmitters selected. 

In these figures the position estimates derived from the differential 
LOP error values are bounded by superimposed LOP segments to indicate the 
range of the LOP values corresponding to the various position points. In 
Figure 2-8 position points are derived from the LOP pair NOR-HAW and TRI-HAW. 
At 10.2 kHz it can be seen that the range of differential error in LOP 
TRI-HAW is approximately symmetric about the zero value (orgin of 
EAST-NORTH coordinates). The points spread into the southwest quadrant 
indicate several error values on the NOR-HAW LOP difference which are- 
positive. The position points are spread along the line which corresponds 
to the azimuth of the TRI-HAW LOP which happens to be in the general direction 
of NOR (29°E). The same general characteristics are true at the other 
frequencies. At 13.6 kHz the range of errors in the NOR-HAW is not as large. 

In Figure 2-9 position points are obtained without using NOR yet the 
larger deviation is along the TRI-HAW LOP direction. This is a geometric 
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Table 2-5. Position Estimate Transformations at Hampton, Va 
for 10.2 kHz, 3 Transmitters (ABD) 


Position 

^Differential LOP 

Error^.. 

Error 

AB AD 

BD 

E 

.0371 -.2297 


X 



E 

-.1701 -.0037 


y 



E 

-.1925 

-.2275 

X 



E 

-.1738 

-.0037 

y 



E 

-.1926 

.0367 

X 



E 

-.1739 

.1685 

y 




B is Trinidad 
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Figure 2-13. Side-by-Side Position Error Considering LOPs NOR-TRI, NOR-HAW 
and HAW-NDK. (a) 10.2 kHz, (b) 11.3 kHz, (c) 13.6 kHz 


problem in that phase errors in the HAW-NDK LOP represent a spatial error 
more than four times as great as the same phase error in TRI-HAW. This 
is because of the distance of Hampton, Virginia, from the HAW-NDK baseline. 
Therefore, even though the differential LOP phase error in HAW-NDK is 
smaller than the TRI-HAW error, the position error range is greater. 

In Figure 2-10 the error in NOR is again predominant with the larger 
position error spread along the TRI-NDK LOP. The reduced error in NOR at 
13.6 kHz is again evident. The extreme position points correspond to a 
slightly larger phase error in the NOR-TRI LOP than was evident in the 
NOR-HAW LOP differential phase difference. 

Considering the least square position error estimate using three 
differential LOP error values, three combinations are illustrated. In 
Figure 2-11 the HAW-NDK geometric dilution of precision error coupled 
with the NOR error is evident at 10.2 and 11.3 kHz. The reduced HAW-NDK 
error and the reduced NOR error at 13.6 kHz are apparent. Very little change 
occurs in Figure 2-12 where NOR-HAW is used instead of NOR-TRI. Figure 2-13 
illustrates a somewhat different situation. The two LOPs involving NOR 
yield correlated position errors forcing the spread along the bisector of the 
angle between the two LOPs. This is also the general direction of the 
gradient of HAW-NDK, compounding the situation. 

Summarizing this presentation, the hourly mean differential error 
lower bound as represented by the SxS data is random within about + .5 km 
based on a position error analysis. During transition periods the error 
can extend to + 1.0 km. With better transmitter coverage in terms of 
signal strength and geometry, this error could be held within +0.5 km. 

2.2 OMEGA Receiver Amplitude-Phase Correlation 

The recorded OMEGA data include amplitude values as measured by the 
TRACOR 599R receiver at the frequency for which this receiver is set. 

Appendix A details receiver amplitude measurement variation as a function 
of input signal strength. These curves have been used in conjunction with 
recorded amplitude values to estimate the received signal strength at each 
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of the receivers for several of the SxS situations. Correlation analysis 
of received phase and amplitude has been used to determine if amplitude 
variations between set-ups might contribute to observed variations in the 
differential phase error mean values. This analysis includes a determina- 
tion of the effect of received signal amplitude on phase measurement 
output for each receiver and provides a comparison of receivers. 

The most significant phase variation in the OMEGA signal is diurnal 
as is the case with amplitude. Of specific interest is short term correla- 
tion between phase and amplitude so that it is necessary to confine this 
analysis to discrete periods of general phase stability within each day 
so as to eliminate the effect of the diurnal. 

A phase-amplitude correlation test was run on several sets of OMEGA 
data taken in a SxS receiver environment. Analysis for each data period 
which may include 3 to 5 days of data is run for the union set of all full 
daytime records and the union set of all full nighttime records within each 
data period. Data from each receiver is analyzed separately. 

Correlation tests consist of scatter plots of phase measurements vs. 
amplitude measurements for each data set as well as a listing of mean and 
standard deviation values for phase and amplitude within the set and a 
calculated measure of correlation 


_ E{ (x - x)(y - y)} 


xy 


a a 
x y 


where x indicates mean of x, a indicates standard deviation of x, E{ • } is 

x 

the expected value operator, and x and y represent the phase and amplitude 

variations with time. The coefficient p is a measure of the linear 

xy 

correlation of the variables x and y, i.e., the nearer in magnitude to 1 
the more linearly correlated. The scatter plot provides a measure of 
linear correlation to the extent that points are collected in a tight 
pattern along a "best fit" straight line. A measure of non-linear cor- 
relation is afforded if the data tend to collect around some non-linear 
curve. 

Amplitude measurements at each 10 second sample time correspond to one 
of four transmitters. These recorded amplitude values are converted to a 
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measure of signal strength ( see Appendix A) . Phase values are transformed 
to remove reference oscillator drift as follows. 

Using the strongest signal which is most stable at the Hampton receiver 
site, the NDK (Sta. D) phase values for each daytime period are averaged 
to determine a daytime mean y_^ for the ith daytime record (i=l, 2, N) 

where N is the number of daytime records in a data period. Using N values 
of y., a set of times t. (t =0, t o =8640, t-=17280, etc.) is defined so that 
each value y^ is displaced in time by 8640 10-second time periods. An 
estimate of reference oscillator drift is then derived according to the 
slope of a least squares line through the points {y^, t^} 



where y is drift estimate in centicycles per 10-second period. 

For the entire data period the reference phase drift is removed from 
the individual transmitter phase data according to 


^NEW^ ^OLD^ 10 Ct " C start^ 


where t gtart ®ay be chosen arbitrarily (all times in seconds). ^qLD^ 
represents the originally recorded phase value at time t and <ji^^(t) 
represents the drift removed value used in the correlation tests. 

Independent drift removal is accomplished for phase measurements 
corresponding to each OMEGA transmitter contained in the data period. 

Drift estimation is made using the NDK data only. 

Once the phase data and corresponding amplitude data are transformed 
each data set is used in a correlation analysis for each transmitter at 
the frequency of the amplitude measurements. Separate analyses are run 
for receivers 1 and 2. 

Table 2-6 summarizes results from analysis of two of the data sets at 
the 13.6 kHz frequency. Data set 56 was obtained in late March 1974 and 
data set 91 during April 1975. It can be noted that no consistent 
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Table 2-6. Results of Phase-Amplitude Correlation Analysis for Two Data Sets at 13.6 kHz 



1) Phase values are In cec relative to local oscillator with drift removed. 

2) Amplitude values are In dB relative to "normal 11 daytime North Dakota signal level 





















correlation between phase and amplitude was observed for either daytime or 
nighttime periods. The TRI amplitude was on the order of 18-20 dB below 
that of NDK during both daytime and nighttime periods of data set 56. The 
TRI-NDK LOP mean value was on the order of -0.50 cec. For data set 91 
there was about 14-16 dB difference in the amplitude of the TRI and NDK 
signals. The mean TRI-NDK LOP phase was near 2 cecs. The scatter plots 
and calculated mean and standard deviation estimates for phase and ampli- 
tude are included in Appendix A. These show very close agreement in the 
general characteristics between receivers for any selected period. 

Based on this analysis, if the correlation between phase and amplitude 
were consistently high for the hourly values, a definite relationship 
between received signal amplitude and phase measurement bias could be 
established. Thus when comparing situations where the average signal 
strength is significantly different, a predictable phase measurement bias 
should appear in the differential phase as an error. This analysis does 
not support any relationship between received signal strength variations 
and the observed period-to-period mean differential error variations. 

2.2.1 Pseudo-lane counting with the OMEGA data .— In analyzing OMEGA 
phase data, either LOP data or individual transmitter station data, it is 
necessary to reinsert lane count information so that statistical calcula- 
tions will be valid. OMEGA phase data as recorded is a measure of phase 
in centicycles (cec) modulo 100 at the frequency of interest. All measure- 
ments are samples at 10-second intervals of received signal phase relative 
to the local oscillator phase (reference). Each relative phase measurement 
is some value in the interval (0,100) at the respective carrier frequency. 
When making a statistical analysis of the data, a transition from, for 
example, 60 cec in one lane to, say, 20 cec in the next lane would appear 
numerically as a -40 cec change instead of a + 60 cec shift. In order to 
preserve the correct direction, magnitude, and average value, it is neces- 
sary to track the phase change as it crosses the lane boundary (100 cec = 

0 cec) . This section describes the procedure used in the statistical 
analysis of data to accomplish lane counting or "pseudo-lane counting" as 
it will be termed since it is simply a procedure for post analysis to 
maintain lane identity in the recorded data relative to the starting lane. 
The starting lane is dependent on supplementary information concerning 
either known position or a position estimate of the receiver. 
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With the first order digital phase-locked loop receiver tracking an 
OMEGA signal no more than a 4 cec phase shift can be successfully tracked 
during a 10-second sampling interval. Any phase changes in the received 
signal which are more rapid will inherently cause a loss of lock condition. 
For the stationary receiver this generally is no problem. In fact, when 
analyzing the recorded receiver output with the receiver locked on the 
incoming signal this limitation can be used to advantage. As phase measure- 
ments are interpreted, a comparison of successive pairs of values can 
readily reveal in which direction the loop is tracking. To illustrate how 
this can be used to determine lane changes, assume that two successive phase 
values are 99 and 1 cec, respectively. If the receiver is in lock, this 
must correspond to a lane boundary crossing which means a phase change from 
99 to 101 cec, i.e., a 4-2 cec change instead of a -98 cec change within the 
same lane. From this reasoning it is possible to define a pseudo lane 
counting algorithm. 

For the beginning relative phase value in a data set for a particular 
transmitter and frequency, an arbitrary lane value is assigned. A lane 
count of zero can be assigned and lane values can be shifted to their true 
value by simply adding a constant after relative shifts have been determined 
within a data set. A lane count register is assigned to the initial setting, 
i.e., LCNT - 0. The initial phase sample is then interpreted as 

(j) 0 = SAMVAL 0 4* LCNT 0 * 100 

where SAMVAL is the sample value of recorded phase in cec and cf> 0 is the 
true value of phase within the defined lane structure. The next successive 
10-second phase sample is then determined using 

<j>! = SAMVAL, 4* LCNT, . * 100 
l i l-l 

where i = 1, 2, to the last sample in the data set. Comparing the ith 

sample with the (i-l)st sample and comparing the difference to some pre- 
determined threshold value 
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1^1 = l+i.i - 4>:| > t 

yields the sample to sample change. If is > T then LCNT^ = LCNT^_^ + 

SGN - <j>p where the SGN (•) function has a value of +1 or -1 depending 

on the sign of A<j> . . If Aiji, < T then LCNT . = LCNT . , . Once the lane value 
i i l i-l 

corresponding to the ith phase value is determined then is replaced by 

<j> . = SAMVAL . + LCNT . * 100 

l li 

where reflects the true phase of the ith sample value in whole cycles 
and fractional cycles (cec) relative to the first point in the data set. 

To determine an appropriate value of the threshold T it is necessary 
to consider not only the receiver characteristics in the ideal situation, 
but also the results of receiver noise on the recorded phase values. As 
was stated, the first order loop can respond, in-track, to no more than a 
4 cec input phase change during a 10-second sampling interval. Digital 
noise internal to the receiver and the effect of editing must also be 
considered when defining the maximum possible recorded phase shift from 
one sample to the next within a lane. 

An automated editing computer program is used to pre-process the 
data. Through analysis of raw data, situations have been observed where 
large instantaneous phase shifts do occur in the recorded output due to 
internal receiver noise. These events are of at least two types. One type 
is generally peculiar to one particular loop on any given occurrence. Note 
that there are twelve loops within the receiver, one for each of four 
transmitters at each of three frequencies. Loop assignments are determined 
by the receiver timing and control logic synchronized to the OMEGA format. 

The loop noise normally is recognizable as an instantaneous phase jump at 
a given sample point and then a phase drift back to. the earlier level 
characteristic of a sudden loss of lock and recovery. Recovery is normally 
accomplished within six minutes. The editing algorithm searches for these 
"spikes" in the data and simply removes six minutes of data from the data 
set beginning with the leading edge discontinuity.. These values are tagged 
as bad data and -not considered for further analysis. Figure 2-14 illustrates 
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(a) 



Figure 2-14. Phase of NDK 10.2 kHz Illustrating the Effect of Spike Editing 
Accomplished with a Computer Algorithm: (a) Phase history with 

spike removed during hour 13, (b) Original phase history 

recorded at Yorktown, Va. 
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an event of this type with a plot of the phase data with a "spike” and a 
plot of the data after editing. The edited data simply appears as a gap 
in the plot. There are situations when such phase spikes have occurred 
in succession so that the edited data may have more than a six minute gap. 
Within a six minute period a uniform input phase change of approximately 
144 cec could conceivably be tracked by the loop (using 4 cec per 10-second 
sample). Thus, if successive phase changes are compared in a lane counting 
algorithm, it is necessary to make the sample-to-sample comparison threshold 
T large enough to allow for within lane changes over at least a six minute 
period. For the stationary receiver the input relative phase should never 
change more than 10-15 cec even with appreciable reference oscillator drift. 

A second type of noise in recorded phase data generally affects all 
loops. The effect appears as if an instantaneous shift in reference oscil- 
lator phase takes place. Sometimes those phase shifts are recovered, but 
more often they are not. The recorded phase suddenly experiences a step 
shift without any apparent loss of track. At some later time the phase may 
suddenly shift back to an earlier level. Figure 2-15 illustrates two situa- 
tions described. One involves a recovered situation and one does not. These 
shifts are normally less than 10 cec but can be up to 50 cec. Note that any 
shift greater than 50 cec in one direction can be interpreted as a shift of 
less than 50 cec in the opposite direction since no lane counting is inher- 
ently a part of the receiver. In any continuous analysis of individual 
transmitter phase data, it is necessary to detect these apparent reference 
phase shifts and correct for them. 

It should also be noted that some of the early phase data was taken 
with the receiver loops configured as second order digital loops . In the 
second order loops, the output phase could conceivably track up to a 29 cec 
phase shift from one 10-second sample to the next. This does not, of course, 
affect expected transitions of input phase but does affect the threshold 
criterion used to detect and edit "spikes" in the data. 

To provide for reasonable capability to track the lane changes in 
relative OMEGA phase when there are periods of missing data through editing, 
it seems reasonable to use a threshold value T - 50 cec. Thus if | A<j> | >50 
the lane count is changed, otherwise the lane count remains the same. Even 
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Figure 2-15, (a) Phase History of TRI 10.2 kHz at Hampton, Va. Illustrating 

the Effect of Apparent Reference Phase Shifts. A shift during 
hour 18 is essentially recovered in hour 20 without a lane 
error. 

(b) Phase Histories of NOR, TRI, and NDK 13.6 kHz at Rich 
Mountain, N.C. Illustrating an Apparent Reference Phase Shift 
During Hour 14 Which is not Recovered. A lane error would be 
introduced in the NDK and NOR phase at the time of this shift. 
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with this choice there are situations which will render lane changes invalid. 
These can be attributed primarily to situations where the reference oscil- 
lator phase changes instantaneously. Further editing is necessary to correct 
for these. Figure 2-16 illustrates a situation in which a valid lane change 
is tracked correctly. Figure 2-15 (b) illustrated a situation where a 
reference phase jump caused an incorrect lane change. 

This pseudo-lane counting algorithm is designed to yield phase data 
which is consistent on an absolute scale within a data set. When calculating 
statistical summaries on recorded phase data and on LOP data, this is neces- 
sary so that calculation of mean values are valid. Even with lane counting 
incorporated it is necessary to visually inspect the data to insure correct- 
ness. It would require extremely sophisticated and complex routines to 
yield better results than the method described and even then absolutely 
error-free lane counting would not be guaranteed. 


r 
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Figure 2-16. Phase History of NOR, TRI, and NDK at Hampton, Va. Illustrating 
Lane Counting in the Data Which is Valid. 
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CHAPTER 3 


ANALYSIS OF DIFFERENT OMEGA ERROR 

Data accumulated over the duration of the ground-based experimental 
program were reduced to hourly mean phase values for further analysis. 

One primary reason for this is simply the magnitude of the total data 
set makes it very difficult to process ten-second samples. Also by 
using data averaging the low-frequency variations in differential OMEGA, 
error can be analyzed essentially independent of any high frequency effects 
which may be receiver dependent. For navigation, particularly in an 
airborne environment, a maximum instantaneous type error description will 
not be available. However, these are very often largely a function of 
the receiver, its time- constant, and the effects of receiver vehicle 
movement. The time constant of the NASA-Langley reciever is by design too 
large (-60 seconds) to gain information concerning instantaneous error 
information relative to the airborne environment. On the other hand, 
general characteristics of the OMEGA error can be well analyzed. The 
essence of the analysis to be presented involves evaluation of individual 
LOP differential phase measurement errors, position estimate errors based 
upon pairs of LOP measurement, i.e.j three station OMEGA data, and position 
estimate errors using least-square methods of combining LOP measurements 
from four station OMEGA data. Generally, mean, rms and CEP (circle of 
error probability) summaries of the various errors are used. In the 
parametic analysis of differential OMEGA receiver separation, range is 
of particular interest. Presentation of much of the analysis is segregated 
by receiver site location for the remote receiver relative to the Hampton 
location. Three frequency results are presented so that comparison is 
made at the various OMEGA carrier frequencies. 

3. 1 Differential Mean 

The differential mean phase errors were calculated from ten second 
differential phase error measurements. Hourly mean values were initially 
calculated. These were averaged over various periods of the day for each 
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data set according to whether the data hour was transition period, full 
daylight, or nighttime. Full data set mean values were also calculated. 

Mean values have been segregated on the basis of LOP station pair and 
frequency. Values obtained from the Tracor phase measurements are treated 
separately from those with the digital receivers. In each of Figures 3-1 
through 3-6 the independent variable is the square-root of separation 
range between Hampton, Virginia, and the respective remote receiver site. 

At each range value, one point is plotted as the mean of each available 
data set at the location. The curves represent connected straight line 
segments between the mean error considering all data sets available at 
each site. 

In Figures 3-1 through 3-6 the leftmost range value in each of 
these plots corresponding to LRC is AR=0 or side-by-side (SxS). The 
spread in these mean values from data set to data set is generally less 
than 5.0 cec at any frequency but appears greater in the full daylight 
condition than for other time periods. For other ranges this time of day 
characteristic is not generally true. It should be noted that only one 
data set was available for six of the sites: FOU, CHA, WAN, MID, BLC, and 
WEE. Therefore, no information relative to variations in the mean is shown. 

For large AR, the RMT location was used for several data periods in Figure 3-1. 
For LOP AB there is a definite negative mean phase error and the magnitude 
correlates with frequency. It can be noted that this error is generally 
larger in magnitude than at closer ranges, although the sign of the mean 
error does fluctuate for closer ranges. This same effect may be observed 
in Figures 3-2 and 3-3 also. With LOP BD (Figure 3-5) the mean values 
are generally smaller, which can be expected with the stronger OMEGA signals. 
Also, the RMT mean values are positive except for 13.6 kHz. This appears 
to be the result of a modal interference shift at 13.6 kHz in North Dakota. 
Similar effects are observed at BDF, NCC and THV, which are all very near 
a critical distance from North Dakota. Similar negative means are observed 
at these sites for LOP BC in Figure 3-4. However, with LOP BC there 
exists more correlation between frequencies. The nighttime shift of 
North Dakota on 13.6 kHz is quite apparent in Figure 3-6(c) on LOP CD. 
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Figure 3-1. Data Set Mean Differential Error vs. Separation 
Range for Norway-Trinidad LOP. 
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Figure 3-2. Data Set Mean Differential Error vs. Separation Range 
for Norway - Hawaii LOP. 
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Figure 3-6. Data Set Mean Differential Error vs. Separation Range 
for Hawaii - N . Dakota LOP. 
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Figure 3-6. Continued., 








Although the data are more limited for the Tracor receiver and the 
variation in the data set means is larger, there is good correlation with 
the DPLL receiver data. For all LOPs the mean error within a 50 mile range 
is almost always less than 3-5 cec, which is a significant result. Beyond 
that range, the inability to predict the mean due to dispersion, local 
effects, and primarily decorrelation error introduces significant variation 
in accuracy. Further discussion will follow relative to navigation position 
errors with differential OMEGA, which probably are more relevant. 

In Figures 3-7 through 3-12 the standard deviation of hourly differential 
mean values for each data set is shown. These represent full period 
statistics and are plotted against differential range square root for each 
LOP and frequency. The solid lines are straight line segment connections 
between the overall mean value of this standard deviation value at each 
site. Two locations have predominantly large variations on all LOPs: 

Ahoskie and Bodie Island. This large variation is apparent at all frequencies 
which would seem to prelude a local effects problem, such as power line 
interference, as a contributing factor. For LOPs involving Norway, the 
variation in mean values for the side-by-side situation is on the order of 
3-5 cec at each frequency, with the 11.3 kHz frequency showing the 
smallest standard deviation values and the least variation. With the 
Trinidad-Hawaii LOP all three frequencies show variations within 2 cec for 
side-by-side with an outlier at 11.3 kHz. For Trinidad-N. Dakota, the 
13.6 kHz shows the least variation while for Hawaii-N. Dakota, the 11.3 kHz 
shows the smallest variation. 

Overall, there is some increase in standard deviation with separation 
range, but it is not very significant. In fact, the Midway location was 
an unusually low standard deviation for all LOPs. There is only limited 
data available for this site however. The Tracor receiver data is not as 
extensive as for the other receiver because only one frequency was tracked 
during each set-up. It is difficult to conclude that the Tracor 
measurements have less mean variation than was exhibited by the digital 
receiver. 


o 


59 



CENT I CYCLES AT PARTICULAR FREQUENCY 




TRACOR 

10.2 



TRACOR 

13.6 



2 o iszii. za C' 2 5 


Norwav-Trinidad LOP Stan- 
dard Deviation of Hourly 
Differential Mean for Each 
Data Set. 


TRACOR 

10.2 


TRACOR 

13.6 



Fig. 3-8. No way- Hawaii LOP Stan- 
dard Deviation of Hourly 
Differential Mean for 
Each Data Set. 
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Fig. 3-11. Trinidad-N. Dakota LOP 
Standard Deviation of 
Hourly Differential Mean 
for Each Data Set. 


Fig. 3-12. Hawaii-N. Dakota LOP 
Standard Deviation of 
Hourly Differential Mean 
for Each Data Set. 


3.2 Differential Average Hourly Error Standard Deviation 

The differential phase error analysis included calculation of an 
hourly standard deviation for each hour at each receiver location. These 
standard deviation values were averaged over various time periods to arrive 
at daytime, nighttime, transition period, and full 24-hour period average 
hourly standard deviations. These values have been organized by LOP and 
frequency, and are presented as a function of the square root of separation 
distance between the respective receiver sites and Hampton, Virginia. 

Plots of these analyses are presented in Figures 3-13 through 3-18. These 
presentations represent the effects of statistical decorrelation error and 
inherent receiver repeatability type error. Individual data points at the 
ranges indicate the average for a particular receiver set-up (data period) , 
while the curves are straight line segments interconnecting the mean of 
each of these data points at the respective separation range. 

There is a general trend for the average hourly standard deviation 
of differential phase error to increase with receiver separation range. 
However, the increase with range is not significant. The linear variation 
of these plots does indicate an exponential relationship between 
differential phase error standard deviation and receiver separation range. 
Local effects are quite evident in these presentations and appear to be 
a dominating effect. The Wallops Island, Virginia location is consistently 
very quiet in that the average hourly error standard deviation is 
generally less than all other locations except side-by-side. The side-by- 
side standard deviation is a good measure of the uncorrelated receiver 
error and is generally less than 1 cec at each frequency with the digital 
receivers. With the Tracor receiver this error is larger on the average, 
generally between 1 and 2 cec. Signal- to-noise ratio is an important 
consideration with the Tracor receiver since the LOPs involving Norway 
(A) generally have larger side-by-side errors. It can be noted that the 
BD LOP errors on Tracor compare more favorably with the digital receiver, 
particularly during the daytime at 13.6 kHz. These are the two stongest 
signals and with the digital receiver there is very little variation in 
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(a) Transition Time (b) Daytime 

Figure 3-13. Data Set Differential Average Hourly Error Standard 

Deviation vs. Separation Range for Norway - Trinidad LOP 
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Figure 3-14. Continued. 
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(a) Transition Time (b) Daytime 

Figure 3-16. Data Set Differential Average Hourly Error Standard^ 

Deviation vs. Separation Range for Trinidad - Hawaii LOP 
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Figure 3-16. Continued. 
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Figure 3-17. Continued. 
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the error with separation range. This indicates good correlation between 
received phase values at the two receivers, even for separation ranges of 
600 n.mi. 

In general, the slope of the Tracor receiver error is greater than 
for the digital receivers. The reason for this is not clear. One con- 
tributing factor is that the set-up to set-up variation (receiver re- 
peatability) is markedly worse with the Tracor receiver. With more 
variation in the individual data points the curve of the means appears 
steeper. 


3.3 Differential Position Error 

Using the transformations discussed in Chapter 2, both two LOP 
(three-station) and three LOP (four-station) differential position error 
has been investigated. It should be noted that with three stations a 
unique position estimate is obtained. With four stations the position 
estimate is not unique. The error is dependent on the LOP choice. For 
four-station data several LOP triplets have been used for some of the 
following analysis and the position error illustrated is the smallest 
in a least-square sense. 

In Figure 3-19 several data sets were combined to yield a set of 
hourly differential position error estimates using the Norway, Trinidad 
and North Dakota phase measurements at 10.2 kHz. The RTI site differen- 
tial data were used corresponding to a receiver separation range of 
258.9 km. The RTI site is in a rural area with the whip antenna mounted 
on the roof of a one story building. The position estimate RMS value 
was 0.60 km with a CEP value of 0.41 km. Figure 3-20 is a histogram of 
the magnitude of the position error values. It may be noted that the 
position errors are relatively evenly distributed directionally and yield 
quite good accuracy at this separation range. 

To compare ordinary OMEGA accuracy (No SWC) with differential ac- 
curacy, Figures 3-21 through 3-26 illustrate RMS and CEP position error 
using stations Norway, Trinidad and North Dakota at each of the carrier 
frequencies and each of the 13 primary receiver sites. These plots are 
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11.3 kHz 



Figure 3-21. OMEGA RMS Position Error (kin) at Various Receiver 
Sites Relative to Hampton, Virginia Site. No SWC 
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Figure 3-22. OMEGA Position CEP (km) at Various Receiver Sites 
Relative to Hampton, Virginia Site. No SWC. 

LOP's NOR-NDK and TRI-NDK. 




13.6 kHz 



Figure 3-24. Differential Position CEP (km) vs. Distance 
Between Receivers (n.mi.). 

LOP's NOR-NDK and TRI-NDK. 




10 . 



S33I3Hni3 


33 






based on RMS and CEP position error over all hours in several data sets 
at each receiver site. The range values are relative to Hampton, Virginia, 
which is at a zero range value. Comparing ordinary to differential, 
the effect of the differential correction is graphically apparent. For 
receiver separation of up to 278 km (150 n.mi.) the differential 
correction offers a significant improvement in position error. Beyond 
that range there was little, if any, noticeable improvement. There 
appear to be some local effects contributing to position error at both 
AHO and BOD, which exist in both the ordinary and differential position 
error plots. The effect of frequency on the error is not clear, although 
the ordinary OMEGA error appears very slightly larger at 10.2 kHz at 
virtually all sites. In Figure 3-27 the overall mean differential posi- 
tion error is shown at each site. The 10.2 kHz mean at BOD is particularly 
large and illustrates further what is probably a local effect problem. 

Note that the RTI location which was addressed in Figures 3-19 and 3-20 
is particularly good. Considering four-station data. Figures 3-28 through 
3-30) compare ordinary (No SWC) and differentially corrected position 
errors using the minimum least squares position error LOP triplet at 
each site. With uncorrected OMEGA there is essentially no advantage 
using LOP triplets over the LOP pairs used in Figures 3-21 and 3-22. 

However, with differential corrections, an advantage does appear. 

3.4 Summary 

Recognizing that any analysis of differential OMEGA accuracy is 
dependent upon the receiver equipment, this analysis does confirm the 
differential concept as a means of providing for usable accuracy. Certainly 
the differential correction can be valid with up to 275 km (150 n.mi.) 
separation of receivers under many conditions. Even at extended separations 
beyond this, good accuracy is possible, as demonstrated with these results. 
This analysis considered only real-time corrections. Other studies have 
indicated that significant time delay in corrections is costly in terms 
of accuracy. Even so, with real-time corrections, position estimate 
accuracies do vary with time of day, and based on these results, are 
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13.6 kHz 
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Figure 3-27. Mean Differential Position Error (km) vs. Distance 
Between Receivers (n.mi.). 

LOP's NOR-NDK and TRI-NDK. 




13.6 kHz 



Figure 3-28. OMEGA RMS Position Error (km) at Various Receiver Sites Relative to Hampton 
Virginia Site Considering LOP Triplets NOR-TRI , N0R-1IAW, HAW-NDK ; NOll-TRI , 
TRI-HAW, HAW-NDK; NOR-HAW, TRI-HAW, HAW-NDK. 
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Values Correspond to Triplet with Minimum Least Squares Position Error. 
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Figure 3-30. Four Station Differential Position Error with Norway-Trinidad 
Trxnidad-llawaii, and Hawaii— N. Dakota LOP's. 


influenced significantly by receiver repeatability. It is very difficult 
to draw conclusions concerning local effects problems and how these affect 
repeatability error. Repeatability error was significant even for the 
side-by-side tests, when the receivers and associated antennas were 
co-located. 


i 


90 



CHAPTER 4 


VLF PROPAGATION 


4.1 General 

A VLF propagation model was developed to identify and better under- 
stand those phase error mechanisms that are related to propagation 
anomalies such as modal interference. Although a complete analysis of 
phase data requires a statistical approach, a deterministic model of 
propagation under average conditions is required to utilize known physics 
to explain important error mechanisms and to provide predictive informa- 
tion to guide the statistical analysis. 

4.2 The Propagation Model: Development of Theory 

Intrinsic to the waveguide theory of VLF propagation is the represen- 
tation of the EM fields as a sum of propagating modes. This mode sum is 
made tractable by introducing a notation to represent the different 
functional characteristics of the solution, as defined by Wait (ref. 10) and 
subsequently extended by Galejs (ref. 11). 

For a Vertical Dipole Source, the mode sum becomes 


E r (h,d) = 


-qldse 


rrrr/4 


h /a A sin (d /a) 


E S 1,5 A 6 G e (h ) G e (h) exp {ik S n d) (4-1) 
q q q q s q r o q ±J 


where h is the distance above the ground surface; A is the excitation 

3 

e e 

efficiency factor of the qth mode; G (h ) and G (h) are the source and 

q s q 

receiver height gain functions respectively; is the freespace wave 
number; d is the distance along the earth T s surface from transmitter to 
receiver; Ids is the dipole current moment of the source; a is the earth ! s 
radius; X is the freespace wavelength. The propagation factor for the qth 
mode can be considered as the sine of the complex angle of incidence at 
the ground, which is commonly used in the literature [refs. 11,12,13). The 
superscript e in (4-1) indicates the vertical polarization of the source 
antenna. A waveguide mode as formulated in (4-1) has both TM and TE 
field components, due to the anisotropy of the ionosphere. However, at 
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VLF frequencies, the field components of one type are usually dominant, 
lending to descriptive terms such as quasi-TM and quasi-TE for the indi- 
vidual modes (ref. 14). To characterize the ground level field patterns of a 
VLF source located on the ground, the excitation factor and propagation 
parameter must be evaluated for each important mode. The height gain 
functions for this situation are defined to be unity. 

An approach based on the work of Budden (ref* 12) represents the wave- 
guide in rectangular coordinates and accounts for the effect of earth 
curvature in the airspace by artificially modifying the permittivity 
with increasing height. The ionosphere is treated, in this approach, using 
a planar geometry. This leads to a "transition" region in which the earth’s 
curvature cannot be neglected but which must also be treated as part of the 
ionosphere. A method by which this transition region can be treated for 
daytime conditions is discussed by Gossard (j:ef. 13). 

Another approach, used in this work, represents the earth ionosphere 
waveguide in cylindrical coordinates so as to take into account the curvature 
of the earth in the direction of propagation. The analysis of the problem 
is divided into two parts. The first describes the fields in the airspace 
below the ionosphere, that region between the ground and the height, tu, 
below which the ionosphere can be completely neglected. The second part 
describes the fields in the ionosphere, that is, above region 1 , up to a 
height, h^, where the continuing variation in the ionosphere is taken to be 
homogeneous. The technique used in the ionosphere region could be extended 
to ground level; but since the airspace solution offers a considerable 
reduction in complexity of analysis, the two region approach is used. 

4.2.1 Fields in the airspace .— In analyzing VLF propagation in the 

airspace region, the TM and TE fields are uncoupled when represented in a 

cylindrical coordinate system, and therefore can be treated independently. 

The TM fields are H , E, , E , and the TE fields are E , H, , H , for propa- 

z 9 r z 9 r 

gation in the $ direction. Figure 4-1 illustrates the coordinate system 
within the earth ionosphere waveguide. 

To solve for the fields in the airspace region, it is convenient to 
characterize the upper and lower boundaries in terms of impedance relation- 
ships at each boundary. A smooth, homogeneous surface can be represented 
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IONOSPHERIC BOUNDARY 



Figure 4-1. Earth-Ionospheric Waveguide Illustrating 
Coordinate System 
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by separate surface Impedances for TM and TE incident waves. It has been 
established that an effective surface impedance for the ground provides a 
good approximation to the actual physical situation at VLF (ref. 10) . The 
TM and TE ground surface impedances are dependent on the angle of incidence 
or propagation parameter. The surface impedances in the limit of grazing 
incidence are 


Zg(TM) 


E z (a) 




iuis 

o 

a+iojs 


Zg(TE) 


-Ej.Ca) 

^ (a) 



-imp /\ iuje 

— r- 2 / VL + r 

J-IWS/ 1 0-1 


O 

i(i)£ 


(4-2) 


which are expressed in terms of ground conductivity a. These impedances are 
numerically close in value to the intrinsic impedance of the media at VLF. 


The conditions necessary for a mode to propagate are that the impe- 
dance condition at the ground is satisfied and that similar conditions be 
met at the interface separating the airspace from the ionosphere. If 
assumptions are made, which lead to an isotropic ionosphere (transverse 
horizontal magnetic field) the TM and TE components are also uncoupled 
in the ionosphere. 

In the isotropic situation a TE and TM ionosphere surface impedance 
can be written in terms of physical ionosphere parameters. TM modes which 
solve the TM boundary conditions can be solved for, independently of TE 
modes. The analytical field solutions which involve Hankel functions of 
large order and argument can be approximated by Airy functions or Hankel 
functions of order 1/3. A modal equation can then be written in terms 
of these Airy functions and the TM boundary conditions. This equation can 
be solved iteratively for the propagation parameters of allowable TM 
modes. The TE modes are found in similar fashion. This method may be 
extended for an anisotropic ionosphere by coupling the TM and TE modal 
equations to satisfy the more complex boundary condition at the ionosphere. 

This condition can be expressed by Budden's four reflection coefficients, 
for example (ref. 12). An alternate approach (re£ ll)-used in this work is to transfer 
the TM and TE impedances from the ground to the ionosphere in vertical 
steps corresponding to thin cylindrical layers. The functions representing 
the radial variation of the fields can then be approximated by exponentials. 
The fields at the boundary of one layer are transformed to the next layer 
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boundary by a TM or TE transmission matrix E or H respectively. For the 
TM case the analytic solution yields: 




-V 

~E r (r )' 
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cue r 
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H z (r) 
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( 4 - 3 ) 


where and are Hankel functions of the first and second kind, and of 
order v. The fields can be approximated over a thin layer using: 
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where K 


1’ K 2 ± 2r ^l ] 


„ 2 v + .25 ; v = S k^a 

K — “■ o ° 

m 0 r ; r = the mean radius of the layer, 

m in 


This result follows by representing the radius as average layer radius 

where it appears as a coefficient in the differential equation of the 

analytic solution. For the nth layer of average radius (r + r , _ ) / 2 , 

n n+1 

A and B can be eliminated to give: 
ee ® 
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Similarly a transmission matrix can be defined for the TE field components. 
These results are dependent on the propagation parameter S. 



The transmission matrices can be cascaded by multiplication to form 
an overall transmission matrix for the TE and TM modes which will give the 
TM and TE impedances at the upper airspace boundary for a given mode, 
since the ground impedance is known. For the TM case: 
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(4-7) 


The analysis in the ionosphere region yields coefficients, a. , that 

ij 

relate the TE and TM impedances and Z^> respectively, at the lower 
ionosphere boundary. 
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a ll + ■ 

a 21 + a 22 Z h 
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a 12 
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a 22 Z e 


(4-8) 


These coefficients can also be expressed uniquely in terms of the four 
Budden reflection coefficients (ref. 11). Given the airspace boundary informa- 
tion, it is then possible to solve iteratively for a value of S which 
provides transformed impedances at the ionosphere that fit the impedance 
relation determined by the ionosphere solution. 

4.2.2 Fields in the anisotropic ionosphere .— A thin layer approach 
is also used in the ionosphere to determine the ionosphere lower boundary 
impedance relation. However, here it is not possible to separate the 
TM and TE field components. Consequently, a transmission matrix is used 

which transforms the four fields E, , H , H. , E across a thin cylindrical 

<p z <p z 

layer. The same coordinate system as that used in the airspace is kept. 

The fields in the ionosphere are expressed in terms of two upgoing and 
two downgoing magneto-ionic modes. The individual wavenumbers are obtained 
from an expression analogous to the Booker quartic. For an assumed time 
variation of the form e Maxwell's equations can be written for the 

ionosphere as : 
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V x E = imyH 


(4-9) 


V x H - -i5e [e] • E 
o 


(4-10) 


where = m/y^e^. As in the airspace, the radial dependence is approximated 
as a complex exponential in a thin layer. The field variation in a layer 
is assumed of the form: 


E = E q exp(-iu)t + ikr + iv<j>) 

Equations (4-9), (4-10) and (4-11) are combined to give: 


(4-11) 


0 = [A] E 


(4-12) 


This is satisfied if the determinant of A vanishes. The determinant of A 
is expressed as a quartic polynomial in k which can be equated to 0 and 
solved to specify the radial wavenumber (k) for each of the magneto-ionic 
modes. 

The field components can then be written in terms of these modes as: 
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(4-13) 


where A^ through A^ are the amplitudes associated with each mode. For 
the purpose of determining the coefficients of C , the amplitudes A^ 
through A^ can be taken to be the magneto-ionic mode components of a 
single spatial field component such as given in (4-11). This formulation will 
not allow a transition to isotropic regions as all four field components 
will not be coupled with the one selected. Therefore, to allow numerical 
transition to an isotropic environment, two ionospheric modes are defined 


97 



as components of E . Evaluation of the elements of C then proceeds from 

z n 

Maxwell's equations and equation (4-12). The transmission matrix T for a 
layer is defined in terms of the field components at radius r and r + Ar, 
where Ar is the layer thickness. 
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The transmission matrix T is derived from (4-13) and (4-14) and can be written 
as 
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where the exponential terms were collected to emphasize the dependence on 

layer thickness and eliminate the need for numerical calculation of exponentials 

ik it 

with large arguments such as e 1 . To express the fields at the bottom of the 
ionosphere in terms of the modes at the top of the varying ionosphere a matrix 
Q is defined following Galejs '(ref. 11) with the relation: 


E, (r.) “ 
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" A, 

9 i 
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V r i> 




A 3 

_ 
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(4-16) 


where it is assumed there are only upgoing waves in the region above h^. 
Equation (4-16) can be rearranged to give the boundary relationship for the 
ionospheric surface impedances as 
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or 


(4-17) 


Z 

e 




Z = 
h 




where through a 22 are expressed in terms of the elements of Q as 


a ll q ll q 43 q 13 q 41 

a 21 = q 21 q 43 " q 23 q 41 


a 12 “ q ll q 33 " q 13 q 31 

(4-18) 

a 22 = q 21 q 23 ~ q 23 q 31 


The results of this model based on a thin layer approximation technique 
in cylindrical coordinates have been compared with the Galejs model. Galejs 
(ref.ll) analyzed the ionosphere in rectangular coordinates and accounted for 
earth’s curvature by decreasing the horizontal propagation parameter, S, with 
increasing height for each layer. This was similar to the modification of 
permittivity with height adopted by Budden and necessarily places a limit on 
maximum layer thickness for the approach to account for curvature in an 
acceptably accurate manner. 

The mode parameters, parameters produced by the method of Galejs and 
the method presented here, are in excellent agreement, with a small deviation 
occurring as the layer thickness is increased. It has been determined that a 
cylindrical layer analysis using exponential approximation for the radial 
cylinder functions gives accuracy comparable to that of the Airy function solu- 
tion for a layer thickness of 10-20 km in the airspace. This indicates that a 
determination of shell thickness in the ionosphere region should be based on 
the resolution in available ionosphere profiles and possible consideration of 
unwanted reflection from layer boundaries. 


99 





4.3 The Propagation Model: Computational Algorithms 

The preceding section provides the theoretical basis for a numerical 
model developed to study the effects of modal interference on VLF naviga- 
tion. The model determines the propagation parameters and the excitation 
factor A which are required by (4-1) to characterize the propagation of each 
mode. The mode parameters can be used to evaluate the extent of phase and 
amplitude interference to the primary mode along a particular propagation 
path. The ionosphere, which is considered as a series of concentric 
cylindrical layers, is taken to be homogenous within an individual layer. 
Within such a layer, the electron density and electron collision frequency 
are required to characterize the ionosphere by a permittivity tensor [e]. 
Profiles of these two ionospheric parameters which are arbitrary or of 
exponential form can be used. Reference profiles evaluated by Deeks from 
experimental data have also been included. The orientation of the earth 
magnetic field with respect to the propagation path is specified to the 
model by the magnetic dip angle and the azimuth angle of the propagation 
path with respect to Magnetic East. The earth magnetic field strength 
along with its orientation complete the terrestrial parameters necessary 
to determine [e] for a given layer. The earth boundary is characterized 
to the model by an equivalent ground conductivity and relative permittivity. 

The propagation parameter S is solved for iteratively for each mode 
in two nested stages which require that particular boundary conditions be 
matched. Since the boundary impedence relationship obtained from the 
ionosphere solution is less sensitive to S than the match below in the 
airspace and is also by far the lengthier in run time, the algorithm is 
written to minimize the number of passes through the ionosphere. An initial 
estimate of S is taken from the startup algorithm to be subsequently 
described and is passed to the ionosphere solution along with the physical 
inputs mentioned previously. The parameters A^ through A.^, which relate 
the TE impedance to the TM impedance at the lower ionosphere boundary by 
(4-17), are the results from this solution. This relationship characterizes 
the upper boundary of the airspace so that a refined value of S that 
satisfies the airspace boundary conditions can be solved for iteratively. 
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This refined estimate is used as a new initial estimate, and the process 
is repeated until S converges within preset specification. 

In the airspace for an estimate of S, the ground impedance can be 
transformed radially for the TM and TE case to give the wave impedances 
at the upper airspace boundary Z e '(S) and Z^' (S) , respectively. These 
impedances can be tested against (4-17) to establish match indexes which can 
be written 


V s) 



+ A 12 Z H 
+ A 22 Z h’ 



(4-19) 


V s) 


A 11 ~ A 2lV 

A 12 " A 22 Z E* 


- Z ' 
H 


(4-20) 


Both (4-19) and (4-20) approach 0 as a match is obtained and can be used 
to find the mode parameters with a standard complex root-finding algorithm 
such as Muller's method. A contour, map of |Mg(S)| is given in Figure 4-2 
showing the first two TM modes. 

In the nearly isotropic case Z and Z„ are only loosely coupled. There- 

E n 

fore, the zero finding routine tends to bypass the roots of one wavetype de- 
pending on the criteria used, unless the initial estimates are good to begin 
with. This problem is avoided by refining the rough TM estimates using 
(4-19) and the rough TE estimates using (4-20). These rough estimates are 
obtained by selecting values of S along the real line corresponding to local 
minima of Mg(S) and 
estimates, respectively. 

The number of cylindrical layers in each region determines the layer 
thickness and the closeness of the approximation to the theoretical model, 
given sufficient computing precision. The number of layers and the thick- 
ness of the varying ionosphere, determined by tu and h^, are minimized for 
run-time economics. The ionosphere calculations represent the greatest 
contribution to run-time, hence the number of layers in the airspace can be 
chosen sufficiently large and fixed. The repeated multiplication of the 
layer transmission matrices, especially in the ionosphere, is subject to a 
buildup of numerical error. However, the steps described above to reduce 
run-time are effective in minimizing error buildup. For a machine with 


M^(S) for quasi TM and quasi TE propagation parameter 
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16 place significance, there is wide latitude allowed in the selection of 
an h for which the solution! ; in the airspacerat VLF.yis. -not dependent. This 

N , .7^ “».* 

region of independence from ranges from approximately 110-140 km. 

Results from this model ’we ; re, compared tO;, those of Galejs over a wide 
range of frequency, ground conductivity and ionosphere profiles. The results 
indicate essentially exact agreement for shell widths of 15 kilometers or 
less. Comparison calculations were based on the various ionosphere profiles 
of Deeks. Agreement of the propagation parameter was to 1 part in 10^ for 
phase velocity and one hundredth of a decibel per kilometer for the attenua- 
tion factor. The excitation factors also agreed to within .01 dB. It 
appears that the two versions converge to identical results in the limit 
of zero shell thickness. The propagation model also gives excellent agree- 
ment with published results from the NELC Waveguide program. Results from a 
typical comparison are shown in Table 4-1. The specific conditions for the 
comparison at 15.567 kHz are those used by Morfitt (ref, 14) in a. similar compari- 
son of propagation models. These conditions correspond to mid-latitude 
propagation under an assumed exponential nighttime ionosphere. The agree- 
ment in phase velocity and attenuation is limited only by the tolerance 
that was used for iteration cut-off in the RTI program. The excitation factor 
is subject to somewhat more variation because it depends on the derivative 
of an ionosphere impedance "match" characteristic with respect to the propa- 
gation parameter. A discussion of this with the researchers at NELC disclosed 
that due to the variation in direction and increment size in taking the numerical 
complex derivative they would expect as much as 1-2 db of variation in the 
excitation factors for runs made with different initial postulated propagation 
parameter values. 

Table 4-1. Comparison of Propagation Models 


Mode Identification RTI LONGWAVE NELC WAVEGUIDE 


TMl 

c/v 

.99636 

.99632 


a db/mm 

.27 

.28 


|A| db 

-8.31 

-7.81 

TM2 

c/v 

1.00703 

1.007128 


a 

1.46 

1.508 


I A | 

.04 

.91 

TE1 

c/v 

.99802 

.99800 


a 

1.44 

1.49 


1 A 1 

-20.46 

-20.28 
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CHAPTER 5 

COMPOSITE OMEGA ANALYSIS 

i ■ » ■ 

Composite OMEGA is a mode of operation whereby phase measurements 
at two or more carrier frequencies are combined to form a composite 
frequency phase. In reference 1 this concept was discussed considering 
composite phase as a linear combination of carrier phase measurements. 

One problem that was addressed involved mapping composite phase values 
to position locations. With no constraints on the carrier phase weight- 
ing this can result in ambiguities. In this chapter weighting coefficient 
constraints are considered and related to work done by Pierce and Baltzer 
(refs. 15, 16, 17, and 18), so that the composite phase may be directly 
related to standard charts. A second section discusses determination 
of combination weights to minimize measured composite phase variations. 

5.1 Composite OMEGA Weighting Coefficients 

Assume that a composite OMEGA phase measurement consists of a linear 
combination of 13.6 kHz and 10.2 kHz phase measurements in units of the 
respective OMEGA frequencies. The composite phase, <J> c , can then be 
represented as 


*1^13.6 + * 2 ^ 10.2 (5 1) 

where (a^, a^) are arbitrary in the range < a^,a 2 < 00 . Phase units of 
the composite measure will be determined by the relationship of a- and 

4 — 1 

a^ so that only when a^ = 1 - j a^ will <j> c be in units of cecs at 10.2 kHz 
(ref. 15). For other (a^, a 2 ) pairs the phase units of will be in units of 
wavelengths at other frequencies. When comparing various composite phase 
measures, e.g., comparing the navigation errors associated with different 
weightings, the units of <j> for any given (a^, a 2 ) pair must be calculated 
because it is necessary to express the various composite phase values in the same 
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phase units. It is, of course, particularly convenient to use cec of 
10.2 kHz since navigation charts are available in. these units. 

For any choice (a^, a^) it is possible to derive a conversion 
factor so that <J>^ is always in units of cec at 10.2 kHz. Let K be such 
that 



(5-2) 


where <f> is as given in (5-1) and 4> f is in units of cec at 10.2 kHz. 
c c 

Then using (5-2) to rewrite (5-1) 



13.6 


+ 


Ka 2 4> 


10.2 


(5-3) 


4 ' 

where Ka 2 = - -j Ka^ is required normalization for <j> c to be in the de- 
sired 10.2 kHz units. This yields 


K = 


[a 2 + 3 V 


-1 


(5-4) 


Thus, the expression for composite phase in 10.2 kHz phase units for 
any choice of (a^, a 2^ weighting coefficients becomes 





> 1 + 13.6 


+ a 2^10 . 2 1 


( 5 - 5 ) 


where Ca^, a 2 ) can be arbitrarily chosen. 

Equation (5-5) can be expanded and rewritten to yield 



13.6 


(m-l)<j> 


10.2 


(5-6) 


which is the composite form used by Pierce (refs. 16 and 17). In 
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(5-6) 


4 

3 a l 


I a i /a 2 


m = 


. 4 

a 2 + 3 a l 


1 + 3 a x /a 2 


(5-7) 


From (5-7) it can be seen that any choice of ( a 2 > a 2 ^ can to a 

particular value of m. In fact, m is only dependent on the ratio of 
^ a l^ a 2^ given in Figure 5-1. This result means that there is 

in effect a many-to-one mapping of the a- , a- plane onto a line a = 

4 1 Z 1 

1 - t a brought about by conversion of composite phase units to cec 
1 

of 10.2 kHz. Actually each line a 1 = 6a with slope 6 and going through 

1 1 4 

the origin is mapped to a single point on the line a^ - 1 - a^ which 

4 

is the intersection of lines a^ = 6a^ with the line a£ = 1 " 3 a l* ^ small 

subset of these lines is illustrated in Figure 5-2. 

When evaluating the composite phase variation as a function of a^ 

and a n all formulations of d) , when transformed to units of cec at 10.2 
2 c 

kHz, will be equivalent for each pair of the set {a^, a^} corresponding 
to a particular 6 value when a^ - Thus, it is only necessary to 

consider variations of composite phase as a function of 6 and the varia- 
tions of composite phase over the entire a^, a^ plane are characterized. 
From Figure 5-1 it was illustrated that for each 6 there is a unique 
value of m. Thus, one can use (5-6) to evaluate the composite phase 
variations as a function of m and consequently will have considered all 
possible pairs of weighting coefficients given in (5-1). 

This observation is illustrated in Figure 5-3. Real data obtained 
at NASA-Langley were used to form composite phase readings. It was de- 
sired to determine what weighting coefficients would yield a minimum 
standard deviation of composite phase when measures of phase from a 
given transmitter at 13.6 kHz and 10.2 kHz were combined linearly. The 
composite phase was calculated as 


*c - 


*13. 6 } - 


B<j> 


10.2 


(5-8) 
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Co®?o sl: - 2 

Weighting 


5-arda-d Deviation vs. H-l 
Phase S<-and<a- ^ Phase. 

f 13.6 kHz and 10- *■ **• 




where phase readings are in cec at each respective frequency. Equation 

3 

(5-8) is equivalent to (5-1) where a^ = — A and a^ - -B. Then, to in- 
sure that <J)^ is in cec units of 10.2 kHz, (5-8) is expressed as in 
(5-5) as 


K ■ '! A *13.6 - B *10.2 } . (5-9) 

Using measurements of 0- - - a ir , and t, 0 from a five day period of 

X J • O , J_vJ • Z XX 

data collected from the Norway station, the standard deviation of 
4> ? (i.e. o’) given by (5-9) was used to calculate and plot a' for a range of 
A,B (0 < A < 6, 0 < B < 6) . The resulting three dimensional plot is 
given in Figure 5-3. The valley of minimum is actually parallel to 
the A-B plane with = 11.7 cec at 10.2 kHz and is along the line 
B ~ 0.6A. The estimated standard deviation from the data for 10.2 kHz 
was 22.9863 cec and for 13.6 kHz was 22.6189 cec. The data estimate 
of T £2 a io 2 a l2 6 was 513 ‘92. It is desired to estimate the value of m 
corresponding to the weighting of (5-6) for which the minimum standard 
deviation of the composite phase is 11.7 cec. 

It can be noted that with B = A-l, (5-9) yields (5-6) with m = A 
and, as has been discussed in this section, the variation of <}/ along 
the line B = A-l will contain all the information about the variation 
of a' c in the A,B weighting coefficient plane. 

Using (5-6) 


<<> 2 = j| “Ms.6 + 


2 - 2 m<m - 1)T 12 a ld.2 c, 13.6 


(5-10) 


With R 12 


T 


12 a iO.2 0 13.6 


I 

a 

c 


( l5 m2 °13.6 + <m " 1)2 ° 2 0.2 ' I 2 


(5-11) 


The mi nimum value of a 1 may be found by differentiating (5-10) with 
respect to m and equating this to zero. This results in 
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_ 2a 


10.2 


- 3/2R 


12 


m o ■ 


9/80 12.6 + 20 W. 


- 3R 


] 


(5-12) 


12 


With a minimum of 11.7 cec the cross-correlation R^ ma Y be expressed 
as a function of m using (5-11) to yield 


3 r_ 

m lc 2 

+ [*5 

my 

2 (11. 7) 2 

(5-13) 

8 m 

-lj 13.6 

3 _ 

m J 10.2 

3 m(m-l) . 


The intersection of (5-12) and (5-13) yields a point (m , R ) which 

l u xz 

allows a minimum a = 11.7 cec. Ideally, the R^ solution should be 
that which was calculated from the data (513.92). However, due to 

T 

errors in estimating a ir . O. _ , , R 10 , and min [a ], there is a slight 

10 . ^ JL j . o 1Z c 

difference. Figure 5-4 is a plot of the two relationships between m^ 
and R given by (5-12) and (5-13). The intersection yields a solution 
of iHq = 2.638. Thus, (5-9) becomes 



3 

4 


(2 - 638) 4>13.6 


1.638(f) 


10.2 


(5-14) 


yielding the composite phase measurement which has minimum variance for 
this data set. For the value of R^ obtained from Figure 5-4, Figure 
5-5 is a plot of vs. m as given by (5-11) which gives the same infor- 
mation as Figure 5-3. The minimum at i 8 2.638 is the global minimum 
for any possible linear combination of <J)^ g and 2 ^ or P art i” 

cular data set. It is of interest to note that this value of m corresponds 
very closely to values observed by Pierce for the Norway signal at Cam- 
bridge, Massachusetts (ref. 17). 

5.2 The Effect of Carrier Phase Correlation on Composite Phase 

It has been shown that -the linear weighting coefficient plane for 
a composite phase measure involving 13.6 kHz and 10.2 kHz is completely 
characterized by the M m-value M weighting of carrier phase values. 
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Figure 5-4. Data Determined R 10 and m for Min [a ']=11.7 cec 





Using the relationship of (5-6) for the composite phase and the expres- 
sion of (5-10) for the variance of this phase it is desired to consider the 
effects of various correlation coefficients. 

First consider = 0, i.e., the carrier phase values are uncorrelated*, 
then from (5-10) 


2 922 ,,.,22 
* 16 m °13.6 + <m - 1) a 10 . 2 • 


(5-15) 


Normalizing by 2 yields 


' 10 . 2 


- _Lm 

16 


13.6 
r 10 . 2 


+ (m-1) 


(5-16) 


Minimizing with respect to m yields 


m = 
o 



1 -1 


+ 1 


as the value of m at which the minimum occurs. 

Reconsidering (5-10) it can be seen that the term involving the correlation 
coefficient R^ then modifies the uncorrelated phase variance so that a minimum 


occurs when R^ = +1 and a maximum occurs when = - 1 • 


Therefore if uncorrelated 


variations in the carrier frequency phase are present they serve to increase 
2 

a c over what would be if variations were perfectly correlated:, making uncorrelated 

noise undesirable. On the other hand if negative correlation exists, uncorrelated 

2 

noise is desirable in that it serves to reduce a . This of course assumes 

c 

m > 1. 


To further present the variation of a c with m consider Figure 5-6 which 
illustrates the minimum a for various values of / °13.6 \ and correlation co- 

efficient R^* Values of R^ between 0.7 and 1.0 are ’given. It can be noted 
that as R ^2 increases for a given value of / CT 13.6 \ the m value, for which the 

' 10 . 2 / 


minimum a occurs, increases 
c 


As the variation in the 13.6 kHz carrier phase 
decreases with respect to the 10.2 kHz carrier phase variation, the 13.6 kHz is 
weighted more heavily in forming the m inimum variance composite phase (smaller 
m) . As correlation between the carrier frequencies increases, the carrier phase 
values are more equally weighted (larger m) in forming the minimum variance composite 
phase. 
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uJard Deviation vs. m for Various Correlation Coefficient 
ferent Carrier Frequency Standard Deviations 


Further analysis of Figure 5-6 shows that when a ^ g and ^ are ver Y 
nearly equal and highly correlated an m value in the range 2 ,5, 4.0 will yield 
minimum This is consistent with results of data analysis presented in the 

previous section. 


5.3 Summary 

This chapter has presented a discussion of composite phase weighting con- 
straints of carrier frequency phase to obtain 10.2 kHz phase units. Furthermore, 
specific data analysis has been used to determine the appropriate weighting 
of 10.2 kHz and 13.6 kHz phase for minimum variance composite phase over a long 
data period. Some additional data analysis of data combined in composite form 
has been carried out by NASA personnel. This primarily involved consideration 
of the 3.4 kHz difference frequency used in lane counting. Generally it was 
determined that the 3.4 kHz is useful for purposes of carrier frequency lane 
determination but is not as good as a carrier frequency for navigational accuracy. 
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CHAPTER 6 


ANALYSIS OF PROPAGATION ANOMALIES 
6.1 General 

This chapter discusses modal interference, a potential source of error 
when using the OMEGA navigation system. This form of error is dependent 
only on propagation path considerations and appears as a phase error in 
the received phase from a particular transmitter. It is not currently 
considered in the published propagation corrections (ref. 19). 

The representation of EM fields at OMEGA frequencies takes the form 
of a sum of propagating modes. One of these modes, the first quasi-TM 
mode, dominates strongly in most situations. Two factors determine the 
importance of a VLF mode at a given distance from the source. Those are 
the attenuation of the mode with distance, a^, and the degree of excitation 
of the mode at the source, A^. In general the higher quasi-TM modes 
attenuate with distance more rapidly than the lower, but can be strongly 
excited, while the quasi-TE modes are poorly excited. Because of the former, 
the use of a station for navigation is not recommended below a minimum 
separation. For propagation west to east at night, the second quasi-TM 
mode can become significant at relatively large transmitter receiver separa- 
tion because of the decreased difference in mode attenuation parameters. 

The purpose of this portion of the study was to verify that modal 
interference had been significant in the data gathered by NASA and analyzed 
by RTI over a three-year period (see ref 1) , and to determine if the effect 
of higher order modes could be predicted accurately enough to improve 
navigation accuracy in the region covered by the experimental program. 

As discussed in Chapter 2, the data base consists of phase and ampli- 
tude measurements for OMEGA signals from four transmitters, NORWAY, TRINIDAD, 
HAWAII, and NORTH DAKOTA. Measurements were made in a differential mode with 
a fixed receiver at Hampton, Virginia, and simultaneously with a movable 
receiver at sites located in the Maryland, Virginia and North Carolina area. 
Ten- second OMEGA measurements were made over 5-10 day set-ups with the 
receiver set-ups repeated at random intervals. Phase, relative to an 
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uncalibrated stable reference, and amplitude were recorded at both receivers 

j 

for the three frequencies and four transmitters involved. 

6.2 Analysis of Data 

Early in the data gathering period anomalous behavior was noticed in 
nighttime LOP phase measurements involving North Dakota transmissions most 
often at 13.6 kHz. This behavior (see ref. 4) involved large phase 
changes during the night in the form of phase shifts between largely 
differing values. Often these changes appeared as level shifts which 
remained for several hours at a time. At other times less coherent large 
scale phase fluctuations were observed. This phenomena appeared to be 
location dependent, occuring at many of the sites but to a much greater 
extent at only a few. This behavior was most clearly observed on the 
Trinidad-North Dakota LOP measurements. Using the OMEGA propagation 
prediction corrections the LOP phase differences can be predicted for the 
Trinidad-North Dakota LOP by: 



where is t ^ ie c ^ art phase difference based on a phase velocity of 

1.0026 times the speed of light, and 4>ppc BD is the difference in propaga- 
tion corrections at the receiver site for the two paths. 

The difference between the hourly mean of measured LOP phase and the 
corresponding prediction was plotted as a function of time of day for large 
segment of the data. From these plots it was apparent that at night a time 
independent mean error was associated with LOP pairs involving North Dakota. 
This was exhibited at most of the receiver locations and was most prevalent 
at 13.6 kHz. Figures 6-1 through 6-9 demonstrate the day to day repeata- 
bility of the nighttime error for the B-D pair at 13.6 kHz. 
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LOP ERROR 2- 93 RIC 5-2 



Figure 6-1. Trinidad - N. Dakota LOP Error Using Propagation Prediction 
Corrections at RMT. Hourly Means at 13.6 kHz. 
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6.3 Predictability of Modal Interference 

The average nighttime LOP phase error was plotted as a function of 
distance to North Carolina with the intention of comparing this plot to a 
theoretical plot of North Dakota second mode phase interference versus 
distance. The B-D nighttime phase error plot is shown in Figure 6-10 for 
13.6 kHz. The data points are plotted as a dot or square corresponding 
to summer or winter. The propagation paths from North Dakota to the 
different receiver sites differ in azimuth from the NDK-Hampton path by 
less than 10°. This small change in path azimuth with respect to the earth's 
magnetic field was not predicted to have a significant effect on the 
North Dakota phase using the VLF propagation algorithm mentioned previously. 
This allows the inferred North Dakota phase to be treated as if taken along 
a single radial from the North Dakota transmitter to facilitate comparison 
with the theoretical plot. The accuracy with which the averaged Trinidad 
signal is corrected at each receiver site affects the agreement possible 
with a theoretical modal interference plot for North Dakota. However, 
Trinidad provides a stable signal in the Hampton area while providing 
propagation path conditions and sufficient distance to minimize the effects 
of modal contamination. The effect of the second TM mode of the Trinidad 
signal in the receiver area is predicted to be less than two cec with a 
maximum site to site variation close to one cec. A plot of the averaged 
nighttime data as described previously, produces the characteristic modal 
interference curve with a slight summer-winter trend noticeable. Important 
also was the steep slope between locations of mean positive and negative 
errors occurring in Figure 6-10 between 1.85 and 1.95 megameters. This 
was originally postulated as the cause of the large time varying phase 
shifts observed at receiver sites at about this distance from North 
Dakota. (ref. 4). The shifts are apparently caused by a "breathing" of the 
propagation parameters resulting in a small shifting of the interference 
pattern. Expected variations in the nighttime ionosphere profiles were 
analyzed to determine if the observations are consistent with the proposed 
explanation. Morfitt (ref. 20) has fit VLF ground stations data to a 
exponential ionosphere VLF model described by 
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-£(z-h-) 

N = N.e 
e 0 

where N e is the electron density at height z , is the reference electron 
density, h^ is the reference height, and B the electron density gradient. 
Parameter changes in hg and B were necessary to fit the data over the course 
of a typical night. Changes of this size, + 20% for B ■ .7 and + 1% for 
h^, have been analyzed with the VLF propagation algorithm, and are quite 
adequate to shift the interference curve so that a receiver location near 
maximum positive phase error could shift to that near maximum negative 
phase error, or vice versa. This can result in a short term fluctuation 
of more than 50 cec at the location in question. Other locations, somewhat 
removed from the large error gradient would not expect to see the large 
fluctuations but would have a significant nighttime mean error as large 
as 20-25 cec, for this particular propagation path. At most of the sites 
studied this mean error appears to be a significant contribution to naviga- 
tion accuracy and yet is not apparent to the OMEGA user as a signal with 
large phase variance. On the other hand, North Dakota provides a strong 
signal and good station geometry in the Atlantic Coast region which in 
the daytime tends to improve navigation accuracy over that obtainable using 
the other receivable stations only. This situation will be further rein- 
forced by the planned phaseout of the Trinidad station. This stems from 
the greater distance and less favorable propagation path associated with 
the Liberia station. 

The quality of the 13.6 kHz North Dakota signal and the general 
consistency of the second mode interference in the Middle Atlantic Coast 
region has led to an attempt to accurately model phase interference from 
higher modes to determine if navigation accuracy could be improved in this 
region. It is expected that the techniques investigated here would be 
applicable to other geographic areas. It should be noted however, that 
propagation path parameters as well as distance, dictate the extent of 
modal interference present so that the West Coast, for instance, can expect 
relatively little high order mode contamination of the North Dakota signal 
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while being at comparable distance. An effort was made to fit the North 
Dakota-Trinidad data to the VLF propagation model mentioned previously. 

The ionosphere electron density profile has a major effect on the degree 
and location of model interference. This profile changes with the time of 
day, the sensor location, and the solar sunspot cycle a representative- 
profile is shown in Figure 6-11. It was determined that the seasonal 
nighttime ionosphere models of Deeks (ref. 21) for sunspot minimum, produced 
too large a summer-winter deviation in the amplitude and location of 
the pattern although a mean ionosphere profile between the two given by 
Deeks fits the data well. Reducing the height of the profile results in 
an improved fit for the summer months, while an increase is necessary for 
the winter months. The changes necessary (+ 1 km) correspond to less than 
half of those indicated from the Deeks profiles. In raising and lowering 
the profiles, the shape was adjusted to approximate the seasonal trend. 

The primary effects of the profile shape on the interference pattern 
result from the height and slope of the main D region shelf at approximately 
85 km above the earth. Varying the height to comply with the Deeks profile 
accounts for the observed direction in the shift of the interference pattern 
in distance. A shift in the shelf slope primarily affects the attenuation 
of the modes, influencing the second mode more strongly than the first. 

This results in a change in amplitude of the interference pattern. 

Height and slope changes with season intermediate to those of Deeks gave 
the best least square fit to the data. The profile shapes in the regions 
and below the main shelf were taken intermediate to the Deeks summer and 
winter figures depending on the height used. The mode parameters and inter- 
ference pattern were not very sensitive to the exact profile shape in these 
regions. 

A theoretical plot of the phase error due to second mode interference 
is shown in Figure 6-12 for a profile of intermediate height between the 
summer and winter extremes and corresponding to nighttime sunspot minimum. 
Figure 6-13 shows a portion of this plot superimposed on hourly average 
experimental nighttime B-D LOP phase error data. It is noted that a 
seasonal profile height variation as great as specified in the Deeks 
profiles produce a seasonal Mode 1 phase velocity variation which is quite 
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Figure 6-11. Deeks Winter Mid-Latitude Nighttime Sunspot 
Minimum Ionosphere Profile 
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large. A variation of this size is not compatible with the published OMEGA 
propagation correction tables which are based in part on experimental data. 
It is interesting to note that the OMEGA skywave corrections obtained for 
1974-1975, the time interval in which most of the data were taken, do not 
show any nighttime seasonal shift in phase velocity. The more recent tables 
do show a slight nighttime shift in the summer months, still much less than 
that which is predicted by the Deeks profiles. 
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CHAPTER 7 


SUMMARY AND CONCLUSIONS 

This ground-based experimental OMEGA program conducted by NASA-Langley 
Research Center has offered a number of contributions to the field. This 
represents probably the largest concerted effort in terms of time duration 
and geographic extent at collecting OMEGA phase data. Although a major goal 
of the program has been the study of the differential concept of operation, 
the data have provided considerable insight into ordinary OMEGA, composite 
OMEGA, local effects problems, receiver repeatability, and propagation 
anomalies. At the outset of the data gathering, the observed modal inter- 
ference phenomenon was actually not expected. The available literature did 
not generally consider modal interference to be a significant problem in 
achieving good navigation accuracies except at positions close to a transmitter. 
To understand observations which contained significant errors caused by modal 
interference - particularly from North Dakota during nighttime periods - an 
unanticipated direction in the research and analysis resulted. It was necessary 
to understand the characteristics of VLF propagation to a much greater extent 
in order to explain why the observed phase perturbations were possible. This 
led to development of a VLF propagation model which retains enough detail so 
that phase variations actually observed could be predicted in terms of propa- 
gation parameters. For the area in which data were collected the model revealed 
a strong likelihood of modal interference with magnitudes of phase shift compara- 
ble to those observed. In the analysis no appreciable confidence in the ability 
to predict when this interference will occur has been obtained. However, to 
alert the navigation community to this problem is considered a significant 
contribution. 

The differential OMEGA results which have been obtained are certainly as 
good or better than for any previous test. The conditions for this program 
were controlled. Phase correlation range was found to be larger than previously 
expected. Good correlation was observed with receiver separation ranges up to 
600 n.mi. In fact, the largest limiting factor appeared as a random perturbation 
which existed even with the receivers co-located . Local effects are certainly 
important. This should not be a major contributor to error in an airborne 
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environment, however. To expect position error to be within the often quoted 
figure of 0.5 n.mi. rms appears to be quite reasonable based on these results - 
even within differential regions extending to 150 n.mi. radius. With good 
receiver-antenna configurations errors on the order of several hundred meters 
are not unreasonable, particularly in daylight conditions. 

Although analysis of the composite mode is not extensive, several general 
conclusions can be made. Difference frequency OMEGA is definitely a viable 
means to lane determination. Better position estimation accuracies can be ob- 
tained with carrier frequency operation with some means of providing propagation 
corrections. In general, various carrier frequency phase perturbations are 
highly correlated with the possible exception of modal type interference which 
usually is more evident at only one carrier frequency during a given period. 

The implications for airborne use of OMEGA are quite good. With inexpen- 
sive programmable digital processors and the simplicity of the OMEGA signal 
format results of this effort show that attainable accuracies are quite 
reasonable particularly for the general aviation navigating community. 

These data that have been collected and catalogued can be of continuing 
used to the OMEGA community. Certainly one obvious use is to serve as a means 
to improve published PPC tables and algorithms for stand alone systems. These 
data represent really a first look at navigation with the 11.3 kHz carrier which 
has received little attention in the past. 
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AMPLITUDE MEASUREMENT IN THE NASA OMEGA RECEIVER 

Amplitude measurement of the received OMEGA signals is accomplished 
using the TRACOR receiver amplitude output signal. At the frequency 
selected on the TRACOR receiver for a particular receiver set-up 
this amplitude output is fed into an A/D converter and scaled so that an 
integer value in the interval (0, 100) is recorded on magnetic tape for each 
received OMEGA transmission. Only four transmitters may be monitored. A 
measurement value represents the amplitude of the respective OMEGA trans- 
mission in dB relative to the normal level of NDK at the Hampton, Va., 
receiver site. There is considerable range in the relative signal strengths 
of stations A, B, C, and D in the differential region used during the 
data gathering phase of the experimental program. Each receiver must 
measure amplitudes which may differ by as much as 50 dB between the level 
of North Dakota and Norway. The amplitude measurement characteristic of 
the receivers is not uniformly linear over this range and therefore must 
be considered when the recorded signal strength data is used to extrapolate 
the actual received signal strength at any given sample time. 

This Appendix presents the receiver transfer characteristics for each 
channel at the 10.2 and 13.6 kHz frequencies used during the data gathering. 

A transformation to convert recorded values of signal strength to received 
signal strength are derived for each channel. The numerical values associated 
with parameters in the transformation are extrapolated from calibration curves 
associated with the amplitude measurement circuit. These calibration data 
were obtained at the end of the data gathering phase and represent the hard- 
ware configuration beginning with data set 56. Prior to the time of data 
set 56, the A/D gain was somewhat larger and the overall amplitude measure- 
ment transfer characteristic was undesirably non-linear for received ampli- 
tudes at the NDK signal level. 

Recovered amplitude measurements will be used to determine if correla- 
tion exists between received amplitude and measured phase in the experimental 
data. 
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A.l Receiver Amplitude Measurement Characteristics 
Each of the TRACOR OMEGA receivers has an input/output transfer character- 
istic which is a function of an operator controlled gain setting on the front 
panel. This control determines the amplitude measurement outputs for a given 
signal level input. Internal to the OMEGA phase measuring and recording 
equipment is an analog-to-digital circuit with a preset non-operator control- 
led gain characteristic. The analog output of the TRACOR circuitry is 
sampled and recorded on tape at 10 sec intervals such that each of four 
channels is sampled once every 4 sampling intervals during periods of data 
recording. 

The output of the A/D converter will be defined as the amplitude measure- 
ment output, an integer in the interval (0, 100), with units of counts which 
are linearly related to a dB value of the signal relative to the "normal" 

NDK signal. The input is signal strength of the particular OMEGA signal. 

The transfer function is the input/output characteristic plotted as output 
counts vs. input signal strength in dB relative to the "normal" NDK signal. 

Figures A-l and A-2 provide transfer function plots of receiver #1 
(fixed base laboratory receiver) channels (A is 1, B is 2, etc.) at 10.2 kHz 
and 13.6 kHz respectively. Figures 3 and 4 provide the same functions for 
receiver //2 (mobile receiver complex) . Each of the figures shows a set of 
curves for each of four gain settings which may be operator set. 

It can be noted that for each receiver the transfer characteristics are 
non-linear. Also the channels within receivers perform differently as 
do the different receivers. 

For each channel, variation of input signal strength is normally within 

a few dB, so that, to calculate what the input was for any given output sample 
value it is adequate to consider each transfer characteristic in a piecewise 
linear fashion. A point on this curve can be associated with a nominal 
signal level of a channel with a given gain setting. The transfer character- 
istic in the neighborhood of this "nominal" point is then approximated by an 
estimate- of the tangent of the curve at that point. Figure '.A- 5 illustrates a 
nominal TRINIDAD level at -25 dB and 11 counts. The linear approximate 
characteristic slope indicates 1.5 dB/count, i.e. a 15 dB variation in input 
signal level yields a 10 count change in the recorded value of signal strength. 
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A. 2 Conversion of Recorded Data to Input Signal Strength 

Let be the recorded value in counts for a particular receiver measure- 
ment of the ith OMEGA transmitter signal strength at the selected frequency. 

Let S^ be the corresponding signal strength of this signal and S q be the 
nominal NDK level or reference. The piecewise linear transfer function of 
the amplitude measurement characteristic is defined in terms of the parameters 
8 , and where (K^, ^i^ defines a point on the transfer character- 
istic in output counts and input dB respectively and is the slope in 
output dB per input count. Then 

[A i - K u ] x + K 2i = 201og S i/S Q (A-l) - 


® i ^2i 

20 ^ A i ” K li J + 20“ 


Table. A-l provides a tabulation of the parameter values for each of the data 
tapes corresponding to side-by-side runs from data set 56 on. These have been 
extrapolated from the calibration data. It should be noted that 
an operator log associated with receiver If 2 (mobile receiver) was used to deter- 
mine the operator controlled gain settings for each data tape. No such log 
exists for receiver If 1 and the transfer characteristic parameter values have 
been estimated using the receiver calibration data and knowledge of the mean 
signal strength from the receiver If 2 data. On the side-by-side runs, the two 
receivers should "see” the same signal so that this should provide a good 
estimate of the parameter values. Included in Table A-l is the frequency at 
which amplitude measurements were made for each data set. 

Figures A- 6 and A- 7 provide plots of average hourly estimated input signal 
strength vs. output recorded signal strength for the Trinidad and N. Dakota 
signals at 13.6 kHz contained on tape #56 for receiver If 1. These represent 
plots of (A-2) using parameter values in Table A-l over two ranges of recorded 
signal strength corresponding to received signal levels of Trinidad and 
N. Dakota at Hampton, Va. 
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TABLE A-l 


AMPLITUDE MEASUREMENT CONVERSION CONSTANTS 


Note: Offset given in counts and dB; slope in dB/count 

To convert RDG(count) - reading in counts - to AMP(dB) - amplitude in dB: 
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TABLE A-l. Continued. 

AMPLITUDE MEASUREMENT CONVERSION CONSTANTS 


Note: Offset given in counts and dB; slope in dB/count 

To convert RDG(count) - reading in counts - to AMP(dB) - amplitude in dB: 
AMP(dB) = SLOPE (d3/count) [RDG(count) - OFFSET (count) 1 + OFFSET (dB) 
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Lgure A-7. Input Relative Signal Strength vs. Output Recorded Amplitude 
for Data Set 1.56 using NDK 13.6 kHz (Ky=0) 


Using the procedure described in this Appendix, received amplitude 
variations in terms of estimated signal strength at the receiver can be 
compared to received phase to determine correlation. With the local 
oscillator phase drift removed. Figures A-8 through A- 13 are scatter plots 
of input signal strength (dB) relative to the mean value versus measured 
phase in cec at the frequency of interest.. Shown in these figures are 
the tabulated mean amplitude with the associated standard deviation, and 
the estimated linear correlation coefficient. Three data sets, 56, 81, 
and 91 are given for both receivers. Data sets 56 and 91 are at 13.6 kHz 
while 81 is at 10.2 kHz. Note that each data set is analyzed in terms of 
daytime periods and nighttime periods for each OMEGA transmitter which 
was monitored. There appears to be no consistent correlation between 
phase variations and amplitude variations for either receiver. 
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(a) Trinidad Daytime, Data Set 56, 13.6 kHz. 

Figure A-8. Input Signal Strength (dB) Relative to Mean vs. Measured Phase in 
° cec with Drift Removed (Base Receiver) . 


156 




(b) Trinidad Nighttime, Data Set 56, 13.6 kHz. 
Figure A-8. (Continued). 
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(c) N. Dakota Daytime, Data Set 56, 13.6 kHz. 
Figure A-8. (Continued). 
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(d) N. Dakota Nighttime, Data Set 56, 13.6 kHz. 


Figure A-8. (Continued). 
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Figure A-9. (Continued) 
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(c) N. Dakota Daytime, Data Set 56, 13.6 kHz. 


Figure A-9. (Continued) 



(d) N. Dakota Nighttime, Data Set 56, 13.6 kHz. 


Figure A-9. (Continued). 
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(a) Norway Nighttime, Data Set 81, 10.2 kHz. 

Figure A-10. Input Signal Strength (dB) Relative to Mean vs. Measured Phase in 
cec with Drift Removed (Base Receiver) . 
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(b) Trinidad Daytime, Data 


Figure A-10. (Continued) 
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(c) Trinidad Nighttime, Data Set 81, 10-2 kHz. 
Figure A- 10. (Continued). 
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(d) N. Dakota Daytime, Data Set 81, 10.2 kHz. 


Figure A-10. (Continued). 
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(e) N. Dakota Nighttime, Data Set 81, 10.2 kHz. 


Figure A-10. (Continued). 
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(a) Norway Nighttime, Data Set 81, 10.2 kHz. 


Figure A- 11. Input Signal Strength (dB) Relative to Mean vs. Measured Phase in 
cec with Drift Removed (Mobile Receiver). 
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(b) Trinidad Daytime, Data Set 81, 10.2 kHz. 


Figure A- 11. (Continued). 
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Figure A-ll. (Continued) 
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, Data Set 81, 10.2 kHz. 
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(d) N. Dakota Daytime, Data Set 81, 10.2 kHz. 
Figure A-ll. (Continued). 


172 





(e) N. Dakota Nighttime, Data Set 81, 10.2 kHz. 


Figure A-ll. (Continued). 
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(a) Trinidad, Data Set 91, 13.6 kHz. 


Figure A-12. Input Signal Strength (dB) Relative to Mean vs. Measured Phase in 
cec with Drift Removed (Base Receiver) . 
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(b) Hawaii, Data Set 91, 13.6 kHz. 
Figure A-12. (Continued). 
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(c) N. Dakota, Data Set 91, 13.6 kHz. 
Figure A-12. (Continued). 
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(a) Trinidad, Data Set 91, 13.6 kHz. 


Figure A-13. Input Signal Strengt-h (d3) Relative to Mean vs. Measured Phase in 
cec with Drift Removed (Mobile Receiver) . 
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(b) Hawaii, Data Set 91, 13.6 kHz. 
Figure A-13. (Continued). 
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Figure A-13. (Continued) 
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APPENDIX B 


DEFINITION OF DAYLIGHT- TRANSITION-NIGHTTIME PERIODS 


For purposes of the data analysis each 24 hour day has been divided 
into four time periods as described in this discussion. The four periods 
are designated as follows: 


DESIGNATOR 

1 

2 

3 

4 


PERIOD 

Sunrise Transition 
Daylight 

Sunset Transition 
Nighttime 


Generally, these time periods are defined in terms of sunrise and sunset times 
at the two transmitters in a given station pair and the two receivers for 
which differential OMEGA analysis is being made (ref. 22). Considering these 
four positions the "sunrise transition period" is the time between the 
earliest and latest sunrise, i. e. , the time between the full darkness 
situation and full daylight situation. The "sunset transition period" is 
correspondingly the time between the earliest and latest sunset. . "Daylight" 
and "Nighttime" are when the four positions are in either simultaneous 
daylight or darkness respectively. 

These periods are defined for analysis based on calculated sunrise and 
sunset times at each receiver and transmitter site. For each 24 hour day of 
data, four times (GMT) are used to define these periods: SRT, start of 

sunrise transition; DAY, start of daytime; SST, start of sunset transition; 
and NIT, start of nighttime. These times are tabulated in Table B-l for each 
period for which data were taken (a data period is normally seven days or less) 
Using this tabulation of values a descriptive variable is defined for each 
hour in the 24 hour day for each day in the data period according to which 
of the daylight condition time periods the hour is to be included in. This 
variable DTIM has a value (1, 2, 3, or 4) for each hour in the data. This 
variable is assigned for each data hour (DAHR) as follows: 
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if SRT<_ DAHR < DAY, then DTIM = 1 

DAY_< DAHR < SST, then DTIM = 2 

SST<_ DAHR < NIT, then DTIM = 3 

otherwise DTIM = 4 

It should be noted that several anomolous situations arise. In some cases 
either DAY=SST or NIT=SRT meaning that there is either no full daylight or 
full nighttime for that particular period. In the extreme there are 
occasions where the latest sunset occurs after the earliest sunrise (e.g., 
Norway - Hawaii at some times of the year) . Here SRT and NIT have been set 
equal to the average of the two times. Similarly it is possible to have the 
earliest sunset before the latest sunrise so that it is necessary to set 
DAY=SST using the average of the two originally calculated values. There 
are no situations where the transition periods are nulled. 
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TABLE B-l. Continued. 


vr 

\n 

WO 

O 

O 

Sf 


O 

O 

O 

o 

© 

a <o 

O 

o 

o 

o 

CN 

z 

i 

o 

o 

© 

o 

CN 

3~ 

VO 

VO 

vC 

sO 

r- 

= 

rH 

—4 

rH 


rH 

rH 

1 O 

o 

O 

O 

CN 



rH 

rH 

rH 

rH 


CN 

CN 

CN 

CN 

<o 


O 

o 

O 

© 

CN 

Ui 

a <n 

cn 

CN 

CN 

CN 

CN 

z 

CN 

CN 

CN 

CN 

CN 

W 

OS CN 


CN 

CN 

CN 

o 

H 

*— » 

*“) 

© 

-4 

rH 

H 

O 

o 

a 

© 

© 


rH 

rH 

f-4 

>H 

rH 

*y 

O 

O 

O 

wn 

NT 


O 

o 

© 

© 

O 

< <o 

CN 

CN 

CN 

CN 

CN 

s 

CN 

CN 

CN 

CN 

CN 

i-t 

ai cn 

N3 

sO 

H3 

VC 

r> 

H 

rH 

—4 

• r-4 

rH 

rH 

rH 

a 

a 

© 

© 

© 


i—4 

—4 

© 

© 

—4 


CN 

CN 

CN 

CN 

CO 


o 

O 

O 

© 

CN 

a co 

o 

0\ 

a> 

O' 

^y 

z 

<N 

rH 

pH 

rH 

rH 

as 

O CN 

-4 

CN 

CN 

CN 

HT 

Z 

rH 

i-4 

rH 

—4 

rH 

rH 

«N 

CN 

CN 

CN 

o 


O 

© 

© 

© 

rH. 

-T 

co 

tO 

fO 

CO 

vr 


o 

© 

© 

© 

o 

2 

< to 

o 


O' 

© 

to 

S3 

CN 

rH 

rH 

rH 

rH 

as 

C Ml 

vO 

vC 

© 

sC 

O 

Z 

| 

—4 

rH 

© 

rH 

rH 

rH 

<o 

co 

CO 

CO 

© 

1 

O 

O 

© 

© 

rH 

*?\ 

o 

© 

o 

© 

CN 


o 

© 

© 

© 

CN 

SC co 

o 

O' 

O' 

O' 

CN 

H 1 

CN 

—4 

rH 

—4 

rH 

1 

as 

O CN 

© 

—4 

rH 

rH 

CN 

z 

H 

—4 

rH 

rH 

rH 

•Hi 

CN 

<N 

CN 

<N 

© 


o 

© 

© 

© 

rH 

1 

cn 






as 






z > 
o a 

z 

a 


Csl 

© 

p-t as 

< 


-3 

ad 

a 

i- ^ 

2 

s 

a 

3 

© 

< rH 


-h. 

— x 

' — . 

— « •» 

© as 

© 

© 

© 

a 

© 

C > 

as 

as 

as 

as 

as 

© © 

-3 

© 

© 

rd 

© 

as 







© 

a 

< 

a 


Id] 

*y 

sr 

o 

O 

© 

,fis 

C 

o 

© 

© 

© 

1 ? 

rH 

© 

rH 

rH 

© 

t-* 







186 



APPENDIX C 


DIFFERENTIAL CHART VALUES 


Table C-l. Calculated Differential LOP Chart Values 
(Values in cec at Frequency Shown) 
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13.6 

11.3 
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86.79 

77.66 

70.21 

75.17 

99.55 

66.05 

88.38 

LRC-FEU 

36.45 

48.61 

40.51 

18.40 

91.21 

9.33 

81.95 

42.59 

68.83 

LRC-FIS 

94.02 

92.03 

93.35 

60.15 

13.54 

77.95 

66.14 

21.51 

84.59 

LRC-AHO 

89.88 

19.85 

33.20 

72.45 

63.27 

36.05 

82.57 

43.42 

2.85 

LRC-WAL 

43.85 

25.14 

4.28 

48.81 

98.42 

98.68 

4.97 

73.28 

94.40 

LRC-RTI 

84.42 

45.90 

4.91 

41.59 

88.79 

23.98 

57.17 

42.89 

.19.07 

LRC-PGO 

63.65 

84.87 

70.73 

88.38 

84.51 

53.75 

24.73 

99.63 

83.03 

LRC-BDF 

34.11 

45.48 

71.23 

43.90 

58.53 

15.44 

9.79 

13.04 

44.20 

LRC-NCC 

81.46 

41.95 

1.62 

7.01 

9.35 

7.78 

25.55 

67.40 

6.16 

LRC-THV 

10.61 

47.48 

89.56 

48.79 

65.05 

54.21 

38.18 

17.57 

64 . 64 

LRC-BIS 

66.07 

88.11 

73.42 

86.58 

48.78 

7.31 

20.51 

60.67 

33,89 

LRC-RMT 

37.52 

16.70 

30.58 

12.55 

83.40 

36.16 

75.03 

66.70 

'5.58 

LRC-FOU 
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99.69 
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25.67 

LRC-WAN 

82.35 

43.13 
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56.01 

13,34 

9.66 
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10,73 

LRC-MI 
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41.26 
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7.00 

LRC-BLC 

85.10 

13.47 

27.89 

46.96 

62.61 

85.51 

61.86 

49.14 

57.62 

LRC-WEE 

94.86 

59.81 

49.84 

57.93 

77.24 

64.37 

63.07 

17.43 

14.52 




NOR-HAW 

A-C 


TRI-HAW 

B-C 

- 

HAW-NDK 

C-D 

10.2 

13. '6 

11.3 

10.2 

13.6 

11.3 

10.2 

13.6 

11.3 

LRC-YKT 

89.62 

86.16 

88.47 

11.5 

82.01 

1.68 

88.03 

84.04 

86.7 

LRC-FEU 

14.19 

85.59 

4.66 

77.74 

36.99 

64.16 

4.2 

5.6 

4.67 

LRC-FIS 

80.15 

40.69 

0.58 

86.5 

48.67 

7.22 

79.63 

72.84 

77.37 

LRC-AHO 

19.61 

59.49 

66.25 

29.73 

39.65 

33.04 

52.83 

3.77 

69.81 

LRC-WAL 

1.73 

2.31 

68.80 

57.89 

77.55 

64.32 

47.07 

96.1 

30.08 

LRC-RTI 

4.0 

38.68 

48.90 

19.59 

92.79 

43.99 

37.57 

50.1 

75.08 

LRC-PGO 

56.7 

42.27 

18.56 

93.05 

57.41 

47.84 

31.67 

42.22 

35.19 

LRC-BDF 

34.65 

12.87 

94.06 

0.54 

67.39 

22.83 

9.24 

45.65 

21.38 

LRC-NCC 

5.04 

73.39 

27.83 

23.58 

31.45 

26.21 

1.96 

35.94 

79.95 

LRC-THV 

48.11 

30.81 

75.68 

37.50 

83.34 

86.12 

0.67 

34.23 

.78.53 

LRC-BIS 

49.95 

66.6 

55.5 

83.87 

78.5 

82.09 

36.62 

.82.17' 

51.8 

LRC-RMT 

16.13 

21.51 

51.26 

78.61 

4.82 

20.68:- 

96.44 

61.88 

84.90 

LRC-FOU 

2.23 

2.98 

2.48 

'83.82 

78.43 

82.02 

85.90 

81.20 

17.67 

LRC-CHA 

30.05 

40.06 

33.38 

48.09 

30.79 

42.32 

15.01 

20.01 

83.35 

LRC-WAN 

16.46 

55.29 

29.40 

34.11 

12.15 

26.79 

75.54 

.72 

83.94 

LRC-MI 

83.54 

78.05 

81.71 

2.70 

36.93 

47.44 

3.60 

71.46 

59.55 

LRC-BLC 

21.91 

29.21 

57.68 

36.81 

15.74 

29.79 

25.05 

33.40 

27.83 

LRC-WEE 

60.47 

13.97 

11.64 

65.62 

54.15 

61.79 

97.46 

63.28 

52.73 
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APPENDIX D 


ESTIMATED AZIMUTH OF LOPS AT HAMPTON (LRC) 

2 2 2 

Consider the hyperbola defined by x - y = a . The slope at point 
^ X l’^l^ ^ e ^ ne< * taking the derivative 

2xdx - 2ydy = 0 

dy _ x 

dx y (x 1 ,y 1 ) . 

Illustrated graphically as in Figure D-l this slope can be found using the 
bisector of the angle between lines drawn from the point of interest to the 
respective foci. 

This provides a method of estimating the slope of the LOPs at any point. 
The azimuth from LRC to each transmitter of a given pair is calculated using 
the CHART program (ref. 1). The mean of these two values is used to get the 
slope of the LOP at LRC. 


N 



lop bc 

A 7 + A Z o 

C B - 180 estimates slope of BC LOP at LRC. 

A BC ’ 2 

The following table provides a tabulation of these slope azimuths. 

A BCD 

South AZ 208.79° 329.40° 102.78° 126.46° 

N AZ From LRC 28.73° 149.40° 282.78° 306.46° 
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Table D-l summarizes slopes of the six LOPs at Hampton, Virginia (LRC) 
involving station A, B(TRI), C, and D. 


TABLE D-l LOP SLOPES AT LRC 


LOP SLOPE @ LRC (AZ) 


A-B 

89.07 

A-C 

155.76 

A-D 

167.70 

B-C 

36.09 

B-D 

47.93 

C-D 

114.62 
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